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ABSTRACT 


A computer simulation model is developed for use in 
analysing the Navy's Fleet Air Defense problem. The model 
provides for employment of monostatic and bistatic search 
radars in a clear or jamming environment. Use of the model 


is demonstrated by analysing a hypothetical problem. 
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i LNTRODUCTION 


The purpose of this thesis is to develop a computer 
Simulation model for use in analysing the Navy's Fleet Air 
Defense problem, and to demonstrate the usefulness of the 
model with a hypothetical air defense scenario. The model 
is unique in that it provides for employment of bistatic as 
well as monostatic Search radars. A bistatic radar is sim- 
ply a radar system in which the radar transmitter and the 
radar receiver are not collocated. This capability was pro- 
vided to allow for the investigation of the possible advan- 
tages of the use of a bistatic radar receiver in view of the 
extensive capability for passive detection of active radars. 

The model, acronym FADS, for Fleet Air Defense Simula- 
tion, is programmed in FORTRAN IV level G, and is designed 
to run on an IBM 360/65 computer meeting the core require- 
ment of 350 kilobytes. 

FADS is structured along the lines of FLOATS IIB, a 
model currently maintained at the Johns Hopkins University 
Applied Physics Laboratory [Ref. 6]. The significant dif- 
ferences from FLOATS IIB which have been incorporated in 
FADS are: bistatic search radar systems, provision for en- 
gagement of bombers by surface-to-air missiles (SAM), the 
addition of a surface-to-surface missile (SSM) capability 
for the defensive forces, a control logic which requires 
the defensive forces to hold fire until fired on, the in= 


clusion of a dispersed defensive force disposition, and the 
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means for’ constructing up to six distinct. target lists for 
the attacking units. 

FADSs isn suitabler tor analysing both large scale problems 
and one-on-one Situations in a clear or jamming environment. 
In general, the model provides for the description of up to 
90 enemy raiders of 10 different classes, 40 defensive units 
of 10 different classes, 15 classes of radars, 10 classes of 
SAMs, 35 time delay distribution functions, and 18 logic 
switches are included which vary the defensive play of the 
simulation. During the computer debugging phase of the 
model construction, FADS was exercised in all of its various 
modes. Included among these were clear and jamming environ- 
ment, Single and multi threats, area and point defense, 
and monostatic and bistatic search radar systems. 

The particular aren lei used to demonstrate the useful- 
ness of the model is an investigation of the relative ef- 
fectiveness of 100 nautical mile and 50 nautical mile 
non-fire control dedicated SAM systems deployed on picket 
type surface units with bistatic search radar capability. 
This type of SAM system employs active or passive homing 
throughout its entire flight. The measure of effectiveness 
for the respective SAM systems is the rate at which enemy 
air-to-surface and surface-to-surface missiles (ASM/SSM) 
arrive at their target area during the ASM/SSM impact phase. 
This is an indication of the level of the enemy attack sur- 
viving the picket defenses and will ionieteients the dimin- 


ishing of the intensity of the air battle which will be 
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waged in the vicinity of the main force. The primary concern 
is to determine what effect the different missile ranges 

will have on breaking up this focus of an enemy ASM/SSM 
atbrack. 

The analysis was conducted with hypothetical or illus- 
trative values for all systems. Specific performance char- 
acteristics for the bistatic radar system and the missile 
System do not exist. Care was ‘taken to ensure that the hy- 
pothetical values were within the range of existing systems . 
and that the proposed systems' parameters were comparable to 
present system capabilities when appropriate. 

Section II of the thesis contains a general description 
of the flow of events in FADS and a discussion of the princi- 
pie assumptions in the model. The specific inpnuts and re- 
sults of the demonstration problem together with an analysis 
Of the results are contained in Section III. . Sections IV 
and V provide a detailed listing of the inputs required and 
the outputs received from FADS. Conclusions are given in 
Section VI. The appendixes contain background SE Onanestalen 
on bistatic radar systems, a verbal description of the simu- 
jation program, glossaries, and a FORTRAN IV computer program 


meSting. 
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II. BACKGROUND AND DISCUSSION 


A. GENERAL DESCRIPTION OF THE SIMULATION MODEL 

FADS is a large scale store-event simulation of a Task 
Force size surface missile engagement with an enemy air 
attack. The model is capable of representing a threat con- 
Sisting of a coordinated missile attack launched from air- 
craft, surface units or submarines. Although the model 
provides for defensive engagement of the enemy missiles and 
launch platforms with either SAM or SSM, there is no deter- 
mination of defensive interceptor attrition of enemy systems 
during the engagement. If interceptors are to be employed, 
they must be played against the launch platforms prior to 
the start of the game and the enemy forces reduced to the 
appropriate level. The Air Battle Analyzer developed by the 
Applied Physics Laboratory of Johns Hopkins University is a 
suitable tool to determine these interceptor attritions 
(Ret , P11). 

The simulation determines friendly force equipment avail- 
abilities at the start of each iteration. These availabil- 
ities are based on a comparison between user inputs of the 
probability that a given equipment is operational and the 
value of a uniform (0,1) random number which is drawn by the 
computer. The next step in the simulation's iterative proc- 
ess is to determine detections and firing events for the 


friendly forces against the enemy launch platforms. These 
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SAM/SSM fire events are stored in an events list in a time 
sequence. The enemy ASM/SSM launch sequence and targeting 
is computed next. Detections and resultant fire events 
against the enemy ASM/SSM's are then determined and entered 
in the appropriate Sequence in the events list. Friendly 
SAM/SSM intercepts and target hits are evaluated followed by 
enemy ASM/SSM hits and ship kill determinations. The kill 
and damage events are time interactive so that offensive and 
defensive force capabilities are decremented as hits occur. 
The simulation then prints statistics for the iteration and 
starts the next iteration by re-entering the flow at the 
point where equipment availabilities are determined. The 
program stops upon completion of the specified number of 


iterations and cumulative statistics are provided. 


B. ASSUMPTIONS 
i Bastatic Radar System 

in the bistatie detection portion of the simulation, 
the principle assumptions made regarding bistatic radar sys- 
tems are: a) the target's bistatic radar cross section is 
equal to its monostatic cross section; b) the target is not 
located along the transmitter-receiver base line; and c) the 
radar propagation losses are the same as those experienced 
in monostatic systems. 

a. . Hadar Cross Section 

nhe approximation for the bistatic radar cross 

section is made on the basis that, in the limit, the bi- 


static cross section is equal to the monostatic cross 
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section except for the case where the scattering angle is 
equal to 180 [Ref.* Lo}: p.' 590]... The’ scattering’ angle us 
the angle between the path from the transmitter to the tar- 
get and the path from the target to the receiver. It is 
equal to 180° when the target is on the base line between 
the transmitter and the receiver. The fact that the mono- 
static and bistatic cross sections are the same in the limit 
does not imply that they are everywhere equal. It does 
imply that ene range of the values will be equal and thus on 
the average allow usage of the monostatic cross section to 
provide representative values for detection ranges during 
the iterative process of the computer simulation. 
be 4 “Target Locat ron 

The assumption concerning the target not being 
on the base line between the transmitter and the receiver is 
necessary Since range ambiguities result when the scattering 
angle is equal to 180°. This assumption is not unreasonable 
Since the typical track may cross the base line, but will 
not originate on it and remain on the base ates A situa- 
tion in which this could occur would be the case when a 
raider while on the base line launches a missile at a target 
which is either the transmitter or the receiver platform. 

c. Propagation Losses 

The Agente tiled concerning the propagation losses 
is appropriate in that similar assumptions are routinely 
made in arriving at detection ranges for monostatic radars. 


A representative propagation loss for the nominal detection 


16 


a 
- 


: - 7 sees © . - a ia _ 7 a | 


; 4 
“~aa said 4 SeStimaewes ohh newt 







Yoresg ody ros 
j 7 wo - ae 
as JE .yeovieor: ctv. df Jomtst wi aret aes 
oan 


wawited Sn¢f eaat sit no wy sata) -ocs eet Pe 
; ; ; . a 
ycHG i 7 ets J ee? oe me 4 Jha Tt oi) =4 


ins / iv y 258 re ) ' é Eg rata 4Ne- 


7 


(visas 7, 


é eve on - 


range is used in conjunction with the signal to noise ratio 
mequired for.a specified probability .of detection .to arrive 
at a minimum reflected power level which must be received 
insorder for .target, detection to oceur. 

An additional assumption concerning the bistatic 
radar system in the simulation model is that the bistatic 
transmitter is always in an operative status. It is not de- 
stroyed by enemy missiles during «he attack phase nor is it 
evaluated as being down or degraded during the equipment 
status evaluation which is conducted at the start of each 
iteration. Individual bistatic receivers have an assigned 
probability of being operational which should also reflect 
the reliability of their associated transmitter. 


. : 
2. Jamming 





Some major assumptions were made in simulating the 
jamming environment in order to keep the problem tractable 
and to remain within reasonable time limits for the computer 
program. The jamming assumptions will tend to drive the 
solution to conservative estimates for the defensive forces. 
The first assumption made is that the jammer and detecting 
radar have a single gain for their main beam. Thus targets 
not in the main beam are considered to be in a unique side- 
lobe which is a specified number of db below the main beam. 
To determine the eal of the jamming power arriving at the 
detecting radar, it was necessary to approximate the posi- 
tion of the jammer relative to the detecting radar. This 


was done by determining the jammer's position at the time a 
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given target was at the midpoint of its track through the 
radar's detection envelope. In the case of monostatic radar 
thes point was the! target's closest point of approach, | for 
bistatic radar the midpoint in the first detection envelope 
was utilized. The jamming power arriving at the detecting 
radar, given the geometry of this midpoint position, was 
then used to calculate the detection range for the target. 
Additionally, in the computation of bomber detections and 
"Shoot" events, no provision is made for jammers which may 
have been shot down prior to the time of detection. This 
assumption is not required in the case of missile detections 
where detection events for missiles take into consideration 
reduction in jamming levels caused by previous destruction 
of jammers. 

The model is limited to four frequency bands for 
frequency separation of individual radars in the jamming en- 
vironment. 

oa. lLarget Tracks 

The principle target track assumptions are concerned 
with minor anomalies from the real world situations. These 
assumptions do not detract excessively from the simulation 
results, but do ease computation and simplify programming. 

A fixed altitude for the raider and missiles, and the lack 
of course changes on the part of raiders are such simplify- 
ing assumptions. In the case where a raider or missile 
track does not change more than 0.1 nautical mile (nm) in 


either the x or y direction, the target is assumed to fly 
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parallel to the non-changing axis. If the target does not 
move more than 0.1 nm in both the x and y directions, the 
target is assumed to be stationary. Additionally, raiders 
who have completed their missile launches continue on their 
preset course rather than exiting the battle area. Although 
these targets will be engaged by the defensive forces if 
higher priority targets are not present, the resulting tar- 
get hits can be readily "backed out" of the simulation re- 
sults by the user if desired. On the other hand, this 
assumption does allow raiders that are providing standoff 
jamming protection for the missiles to continue to close the 
defensive forces and thereby increase the effectiveness of 
their jamming. 

Other iesSs peneral assumptions are covered as they 
occur in the details of the subroutine descriptions in Ap- 


pendix B. 
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III. DEMONSTRATION PROBLEM 


This section presents a demonstration of the use of FADS 
in analysing a hypothetical problem involving bistatic radar 
systems and surface-to-air missiles (SAM). The problem is to 
select the Te ssiiith range required for the SAM system to be 
employed by a small surface combatant (SSC). The proposed 
SAM system is to protect the major units of a naval task 
force involved in an engagement with enemy missile systems. 
The SSC will have a bistatic search radar and the SAMs will 
be active homing missiles requiring no fire control director. 
As such, the SSC will emit no radiation. The SAM character- 
istics and performance capabilities have been specified with 
the exception of the maximum range. Two range capabilities 
will be considered, a 100 nautical mile missile and a 50 
nautical mile missile. The problem thus becomes that of de- 
termining the relative effectiveness of the SSC when equipped 
with the different weniee missiles. As discussed in the in- 
troduction, the measure of effectiveness used to evaluate 
the missiles is the number of enemy missiles penetrating to 
the force center per unit of time. That is, the number of 
missiles arriving at the force center per minute during the 
time period in which missiles are impacting at the force 


center. 
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A. THE SCENARIO 

The BLUE force consists of one High Value Target (HVT) 
and eight small escorts. ~The BLUE force is maneuvering 
within range of enemy air forces. The HVT is located at 
force center and the pickets are equally spaced on a circle 
one hundred miles from oe HVT. The picket force is silent 
except for one early warning aircraft which carries the 
transmitter for the bistatic radars. For purposes of this 
analysis, the aircraft is assumed to be at force center and 
at 20,000 feet. 

The RED force has detected the radar transmissions of 
the aircraft and the HVT and have launched an attack of 
twenty-one aircraft accompanied by three jammers. The RED 
aircraft wili attack simultaneously in three sections, 
separated by forty degrees of bearing. The RED air force 
has practiced this maneuver to the point that the maximum 
deviation from simultaneous launch will be fifteen minutes. 
Each aircraft will take a maximum of three minutes to launch 
its missiles. . The three sections will approach their target 
from 320, 000, and 040 degrees true. The jammers will carry 
equipment to cover all radar bands used by the BLUE force. 
Each section will consist of one heavy and six light attack 
aimneratt. 


The RED air force will attempt to penetrate the BLUE de- 


mense and sink the HVT. The HVT is considered to be unarmed. 


im order to target a BLUE unit in this dispersed formation, 


a RED unit will have to have positive identification of the 


21 











» ly eee 
; al 


' 
i= Te, Ag? y a e% 6 ' Fr 
pel * ok © 
: Ric % h 
| 3 : 5 fr ig 
A. pe 









if +; ° ~y 
ay ’ ~ 
ra atet 
ie : t 
j aad 1 
7 

7 © a %>* ; 
7 is? 3 - 
| Heal 


mnit as hostile. ,An ECM. fix, or visual ID will be adequate. 
As the BLUE pickets are not radiating, no ECM fix may be ob- 
tained. As the pickets are assumed to be smaller than a 
destroyer escort, and the RED aircraft will fly at a minimum 
altitude of 33,000 feet, no positive visual identification 
will be allowed. 

The radio and radar silence of the pickets leaves them 
‘without the ability to coordinate their defense. No BLUE 
intercept will be allowed if the intercept point is within 
two miles of another BLUE unit's assigned position. The 


BLUE units will use a shoot-look-shoot doctrine. 


B. THE INPUTS 
- Lhe Units 
a. The BLUE Units 

the BLUE unats consist of the HVT, eight piekets’, 
and one early warning aircraft. Positions are as indicated 
ima the scenario. | 

(2) NVI. ‘For the purpose of this analysis the 
HVT is assumed to be unarmed. 

(2) Pickets. The pickets will employ a type 4, 
bistatic, band 4 radar with an antenna height of 50 feet. 
They will carry thirty missiles in canisters. The probabil- 
my OL a correct target identification will be .95. 

(3) Early Warning Aircraft. The BLUE early 
warning aircraft will be unarmed. , It will carry a type 4 
radar and serve as the transmitter for the pickets' bistatic 


receiver. 
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b. The RED Units 


The three classes of RED units are heavy attack, 


ich attack and 


(1) 
erait ‘will ‘carry 


acteristics. 


(2) 


eratt will carry 


jammers. 
Heavy Attack) 1°The* RED “heavy ‘attack air- 


a type 5 radar and have the following char- 


WEAPONS : 1 type I ASM 
4 type II ASM 
SPEED: 415 Knots 
JAMMERS: none 
RADAR CROSS SECTION: 40 square meters 
ALTITUDE: 35,000 feet 
faghtvattack NO Thes RED Light attack air= 


a type 5 radar and have the following char- 


@eteristics. 

WEAPONS : 2 type II ASM 

SPEED: 420 Knots 

JAMMERS: none 

RADAR CROSS SECTION: 13 square meters 

ALTITUDE: 33,000 feet 

(3) Jammers. The RED jammer will carry no 

Badar. It will) have the’ following) Characteristics. 


WEAPONS : none 
SPEED: — 417 Knots 


JAMMERS : bands 1-4 
250 watts/MHZ/band 
359° beam width 
-20 db side-lobe ratio 
jamming target, HVT 
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RADAR CROSS SECTION: 32 square meters 
ALTITUDE: 40,000 feet 
2. The Radars 
The two radars which the simulation will use to play 
the simulation are the class 4 and class 5 radars. Class 4 
radars are found on all pickets and class 5 on all raiders. 
a; " Class “4-Radar 
LYPE: “Bistat ic 
Be 120 nmin? 
Oe 30 am/fm? 
SIDE LOBE RATIO: -20 db 
MAXIMUM INSTRUMENTED RANGE: 200 miles 


RADAR FREQUENCY BAND: 4 


TRANSMITTER UNIT: aircraft overhead HVT 

TRANSMITTER ANTENNA HEIGHT: 20,000 feet 
be Class S Radar 

tYPE: Monostatic 

B: 100 nm /m? 

Q: 50 nm /m2 

SIDE LOBE RATIO: -20 db 

MAXIMUM INSTRUMENTED RANGE: 200 miles 

RADAR FREQUENCY BAND: 3 


BEAM WIDTH: § 3 degrees 
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3. The Reaction Time Distributions 
The reaction time distributions for target process- 


ing by the picket are displayed below. 
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CUMULATIVE PROBABILITY 


Coe 
* 0 | ee 2 a a ae a ae 
O. O5 1.0 iS) 2.0 
DETECTABILITY TO DETECTION DELAY IN MINUTES 
Figure 1 
| Ogeechee eS 





CUMULATIVE PROBABILITY 


O. 0.5 8, 1.5 2:0 
DETECTION TO ID AND DESIGNATION TO WEAPONS EQUIPMENT 


Figure 2 
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CUMULATIVE PROBABILITY 


CUMULATIVE PROBABILITY 





ASSIGNMENT TO FIRING CONSOLE IN MINUTES 
Figure 3 
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FIRING DELAY IN MINUTES 
Figure 4 





4. The Weapons 
a. The BLUE Weapons 


There are two eee of BLUE weapons. Both are 
SAMs and are denoted as SAM LR (long range) and SAM SR 
(short range). 

Ci) eau ER Stns wiki be the LOO mile version 
of the BLUE missile. It will have two sets of intercept 
parameters. The first is for a large, slow, high flying tar- 


get. The second is for a small, fast, medium altitude target. 





Figure 5 


SAM Intercept Contour for High Altitude Engagement 


The input value for the SAM LR class are: 


A = 100 

B = 100 

Ne ne. 0 

SLOPE = -20 

SEEK = -5 

SPEED = 1400 Knots 
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PROBABILITY KILL FOR SINGLE SHOT (PKSS) 
= .88 


PKSS (CROSSING) = 
SYSTEM RELIABILITY 
LAUNCH RELIABILITY 


MAXIMUM ALTITUDE 


MINIMUM RANGE 


The intercept parameters A, B, N, SLOPE, 


SEEK are defined in the description of subroutine TIMMY in 


Appendix B. 


5 (2) 


= 909 


"oo 


59), 000 feet 


= @25 males 





Figure 6 


SAM Intercept Contour for the Medium Altitude Engagement 


A = 100 
B = 90 

N = 2.8 
SLOPE = 


All other parameters are the same as 


-1 


high altitude engagement 
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(2) SAM SR. This will be the 50 mile version 
of the BLUE missile. The intercept contours will be the 
same as those for the SAM LR except the down range and cross 
range values are half the SAM LR values. 

b. The RED Weapons 
The RED aircraft carry two types of ASMs. 
Gi) BASMer oe 4. 
| SPEED: 750 Knots 
ALTITUDE: 20,000 feet 
LAUNCH RANGE: 120 miles 
LAUNCH RELIABILITY: 0.85 
INFLIGHT RELIABILETY. 0.80 
RADAR CROSS SECTION: 1.79 square meters 
(2) ASM Type II. 
SPEED: 900 Knots 
ALTITUDE; 15,000 Sfeet 
LAUNCH RANGE: 110 miles 
LAUNCH RELIABILITY: 0.85 
INFLIGHT RELIABILITY: 0.9 


RADAR CROSS SECTION: 1.01 square meters 
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C. RESULTS 


A summary of the results for 100 iterations is presented 


below. 
1. ASM Results. SAM SR SAM LR 
Available 5100 5100 
Destroyed with launch platform WI 2 C BiG) 246 16 4.84) 
Launched 4928 (96.64) 4854 (95.22) 
Failed at launch 724 C472) FI0NGS 47) 
Failed in flight L29) ( “8\.s) 408) (48.367) 
Destroyed by SAMs 0 0 
Impact on HVT 3775 (76.64) 3686 (75.97) 


Dis SAM Results. 





Available 24000 24000 
Launched _ 3401 6440 
Kills Zion G625274) 250S8 Gie2 
ASM 0 0 
Raiders ZENG 2398 
Overkills 605° G72 87) © 2694 (41.87) 
Failures 680 (20:04) 1348 (20.97) 


ag ASM Impact Intervals (minutes). 


Mean 1.468 2, 432 
Variance 3 0.033 0.037 
Skewness -0.060 0.070 
Kurtosis -0.620 -0.532 
Minimum Ibe (Oj 0.966 
Maximum 1.849 1.862 
4. ASM Impacts per Minute. 
Mean 261 0": 26.192 
Variance W4565i 7s SAA 
Skewness -0.584 0.453 
Kurtosis -0.006 -0.136 
Minimum 18.422 be el 8.21.29 
Maximum 31/2561 37.246 


Dy Histograms of the ASM impacts per minute are shown in Figures 7 
and 8. 
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Histogram of the ASM Impacts per Minute, SAM SR 


Figure 7 
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Histogram of the ASM Impacts per Minute, SAM LR 


Figure 8 
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D. ANALYSIS 
The ahalysis of the results consists of three parts. 
“The first part contains a parametric evaluation of the ASM 
impacts per minute for the two SAM systems. The second 
Part provides a nonparametric evaluation. Part three pro- 
-vides general comments concerning anomalies in the results. 
f Pasamet cic Evaluation 

The purpose of the parametric evaluation is to test 
whether the mean ASM impacts per minute against the SAM SR 
(us,9) is Significantly greater than the mean ASM impacts per 
mamute against the SAM LR (u,5))). With its shorter range 
capability, the SAM SR system should attain a lower number 
of kills than the SAM LR system. Consequently the number of 
ASM impacts on the HVT that occur during the defense by the 
SAM SR system should be greater than the number of impacts 
during the SAM LR defense. 

a. Null Hypothesis 

Hy: Uso < Uji99; SAM SR effectiveness is equal to 

or greater than the SAM LR system's as indicated by the mean 
ASM impacts per minute. 

b. Alternative Hypothesis 

je Uso) 2 Uapng SAMUSR is Less effective than 

SAM LR and therefore has a higher mean ASM impacts per min- 
ute. 

C.> Significance Level 


AVsigniticance Level of 0.105 was used in testing 
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ad. Statistical Test 
The null hypothesis was tested against the al- 
ternative hypothesis using a t test. The t statistic for 


testing equal means is: 


2830 =P 255 0) 
Ee = 
pBsotNioo pbsoSso*MiaoSioo) 
NsoNioo0 nis oF oo= 
where 
<p sample mean for SAM SR. (26. 10) 


X,00 = Sample mean for SAM LR. (26.19) 
s = sample variance for SAM SR. (14.651) 
Stoo = Sample variance for SAM LR. (17.541) 


Nso= Ni00= number of samples for SAM SR and SAM LR. (100) 


The null hypothesis will be rejected if the computed test 
Statistic t is greater than or equal .to the table value for 
the 0.95 point of t distribution with 198 degrees of freedom. 


the table value is: 


e. Decision 
The computed test statistic is t = -0.1578. 
Sance the computed t ratio is less than the tabulated value, 
the null hypothesis cannot be rejected at the 0.05 level. 
i. Conclusion 
The effectiveness of the SAM SR system was equal 
to or greater than the effectiveness of the SAM LR system at 


the 0.05 level of significance. 
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g. Justification for Use of the t Test 

In order to use the t test for difference in 
means, three basic conditions ‘had to be met. These were: 

i) the sample means were independent; ii) the underlying 
distributions from which the samples were drawn was normal; 
and iii) the underlying distributions had equal variances. 

(1) Independence. The independence of the 
sample means was a reasonable assumption since the individ- 
ual samples were obtained from different simulation runs 
using different random number streams. 

(2) Normality. The normality of the sample 
mean values could be asserted on the basis of the sample 
size (100) and the central limit theorem or by demonstrating 
that the means were computed from samples which came from a 
normal distribution. The underlying normal distribution can 
be shown on ane basis that both samples. were unimodal and 
Symmetric, and that the plots of their order statistics 
against the inverse normal probability scale was linear. 

The unimodality of the samples is seen in the histograms in 
figures 7 and 8. Symmetry is evident in both histograms and 
Supported by the fact that the kurtosis of the sample values 
#Ss €lose to zero. The order statistic plots, which are 
clearly linear, are shown in figures 9 and 10. 

(3) Equal Variance. The relationship between 
the variance of the population from which the samples were 
Oeawn 1S tested using the F ratio. This test is applicable 


Sineeé the populations are normally distributed. 
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Inverse Normal vs SAM SR Order Statistics 
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Inverse Normal vs SAM LE Order Statistics 


Figure 10 
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Ca) a sNull Hypothesis, Ht: o%. =.¢o7 the 


0 50 1700) 


population variance (o{,) of the SAM SR system is equal to 


the population variance (co? ._) of the SAM LR system. 


100 


()ieriivernative Hypethests. HH: 62, 7 9, . 05 


the population variances are not equal. 

(c) Significance Level. A significance 
level of 0.05 was used in testing H). 

(d) Statistical Test. H, was tested 


against H, using the ratio of the unbiased estimators of the 


two population variances. 


2 
pe SEO doe! _ 1.1973 


ee 14.651 





The null hypothesis will be rejected if the computed ration 


Pero greater than the tabulated value of Le > ssi aes 
ay 5) ) 


The tabulated values are: 


or less than Lene yo be 


ry 
| 


= 1.491 


ry 
| 


= 0.684 


(e) Decision. Since the computed ratio F 
is not greater than or less than the tabulated values for 
the F distribution, the null hypothesis cannot be rejected. 

(f) Conclusion. The samples were drawn 
from populations which had equal variances. | 

2. Nonparametric Evaluation 
The purpose of the nonparametric evaluation is to 
test whether the sample values for the SAM SR system and 


the SAM LR system are from the same population. 
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a. Null Hypothesis 
H,: There is no difference in the number of ASM 
impacts per minute when defending with the SAM LR system or 
with the SAM SR system. 
b. Alternative Hypothesis 
H,: There is a difference. 
c. Significance Level 
A significance level of 0.05 was used in testing 
Ho. 
Oe olatast real Test 
The Kolmogorov-Smirnov (K-S) two-sample test 


was used to test the null hypothesis. The K-S test requires 


the use of independent samples. Independence was shown in 
D.1.2.(1) abcve. The K-S test is conducted using the sta- 
tiUst Le 


= Pr = 
K max ni) Ga) 


where FOX) is the empirical cumulative distribution func- 
tion for the first set of sample values and Gor) is the 
empirical cumulative distribution function for the second. 
The null hypothesis will be rejected if there is 
a large enough deviation between the two sample cumulative 
Geseributions. ~ The critical values: of D for the level of 
Significance of 0.05 and 100 observations in each sample is 


OPt9233. 
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e. Decision 
The computed value of K was determined to be 
0.0800. Since the computed value of K is less than the 
tabulated critical value, the null hypothesis cannot be re- 
jected at the 0.05 level. A plot of the empirical cumula-. 
tive distribution functions for the two sets of sample 
values iS provided in figure 11. This plot clearly shows 
the close relationship between the two underlying popula- 
tions. 
fi. .Conclusivon 
There was no difference in the defensive capa- 
bilities of systems in the environment in which they were 
examined. 
3. General Comments 
The principal anomalies present in the data include 
the following: 
a) a lack of difference in the SAM LR and SAM SR 
effectiveness; 
b) no SAM kills of ASMs; 
c) a greater number of overkills for the SAM LR; 
d) the number of SAM failures; and 
e) the short ASM impact interval. 
These anomalies do not significantly detract from the con- 
clusions presented in the parametric and nonparametric eval- 


uations. 
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a. Lack of Difference in Effectiveness 
The most probable reason for the lack of differ- 
ence in the effectiveness of the two SAM systems was the 
jamming environment. The reduced detection ranges experi- 
enced in jamming will tend to preclude the use of SAMs at 
the longer ranges available to the SAM LR system. 
b. No SAM Kills of ASMs 
The lack of SAM kills of ASMs was primarily due 
to the geometry of the ASM launch positions and the program 
logic concerning SAM targeting. The program logic requires 
that the closest target be considered as most threatening 
and therefore engaged first. In the problem, the ASM launch 
positions were less than 10 miles from their point of clos- 
Sst approach to the SAM ship positions. Consequently, the 
ASMs with their greater speed were closer than the raiders 
to the SAM ships for only the first 30 seconds after their 
launch. After that period, their speed caused them to be at 
a greater range from the SAM ship. Thus the SAM ships had a 
30 second period in which to detect, acquire, aaa designate 
an ASM as the most threatening target. The reaction time 
distributions for the SAM ships precluded such a short re- 
action time in most cases. 
eG. Overkills 
| The Hees number of overkills for the 100 mile 
SAM was due to the longer flight time of the SAM. This 
coupled with a lack of coordination resulted in a greater 


Saportunity for another unit to obtain a kill on the target 
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while a second SAM was in flight. Targets that had been 
destroyed prior to the scheduled SAM launch time were not 
engaged. 
d. SAM Failures 
The number of SAM failures includes SAMs which 
were fired on the ASMs that impacted on the HVT while the 
SAM was in flight. These attempted engagements occurred 
since the SAM units were not given the capability to deter- 
mine the ASM impact point. 
e. Short ASM Impact Interval 
The tightness of the ASM impact interval may 
appear to be pessimistic in relation to a real world threat. 


However, this should not detract from the comparison of the 
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IV. MODEL INPUTS 


The inputs required to use FADS fall into three groups. 
These are the X vector, the target list assignments, and the 
target lists themselves. The amount of actual input re- 


quired depends upon the size of the problem to be analyzed. 


A. THE X-VECTOR 

All X vector inputs are real values and are read into 
the X vector in the following way. Each input card must 
contain six numbers. The first number is the integer address 
rm the Xevector. “This tnteger should be placed in columns 
1-10, right justified. The remaining five numbers are the 
values to be inserted in the five X vector positions begin- 
ning with the position indicated in columns 1-10 of the input 
card. The format for each card with X vector data is 
CeLOPoFLOUS) Uo (li there are no’ data tocbe read into a posi- 
tion, no input is required, unless the position falls within 
the scope of another input card. The end of the X vector 
impurs iscsignaled ‘by the ‘inclusion of a -zero card. This 
card contains a zero in column ‘10 and zeros in the five data 
positions. A brief overview of the organization of the in- 
put array is provided below. Details concerning the indi- 
Vidual sub-sections are provided at the beginning of each 
sub-section. Specifics dealing with the input values are 
discussed in the description of the particular input value 


ana as footnotes to the individual sections. 
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The followingias the organization of the input X vector: 


general 


KCL) tHosX¢L00),, Control. switches 
XCLO1 ) to) XC 500), Ship Classes 
X(501) to X(900), SAM Classes 


X(901) to X(1200), Radar Classes 


and program constants 


MCLZOL) to X(1400)., Ship Positions 


X(1401) to X(1800), Delay Functions 


X(1801) to X(2000), Raid Classes 


X(2001) to X(2200), Raid Weapon Classes 


(2201) to X(2740), Initial Raider Positions. 


Control Switches 


These inputs include logic switches for varying the 


play of the game. They also provide the means for 


constrolling the output in the form of debug and array 


printout 
lated to 
sections 


portions 


switches. Several of the control switches are re- 


inputs for specific systems which appear in later 


and must be consistent with later usage. Other 


of this anput Section are wsed to store values for 


generating statistics and for internal use by the model. 


These portions should not have inputs since their value will 


be set as required by the model. 


ko 


Inputs 


Case number 


Blue force equipment availability switch: 


O = 100% availability 
ab 


= use input availabilities 
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12 


13 


Inputs 


Blue force radar detection range printout in 
Event 2: 

= no prant 

iy = print 


Fleet coordinationt of Biue forces: 
O = none 


1 coordinated defense 

Blue shrp toubofiiaction after N hats: 
OQ = yes 
1 = no 


Qieais) a Shap iclass) aapuit ) 


Number of iterations. (Must be at least 1. If 
X(85) used, must be at least 10 and less than 
HOA4)) 


Horizon/Sector Scan for Blue fire control radar 
system enabled: 

0 no 

ib yes if equipment capable 


tt oll 


Blue SAM Sector enable: 
+ 


—_ Gone use sectors 
waa v wow wuUwy Vaew 


nN 
1 = sectors used 


Pranto@ontror: : 
O = none : 
1 = selected output tables and intermediate 
debug printouts. Set to O by the model 
aLbpen £inrsth iteration 


Detectability to detection delay: 
0 = use input distribution function 
OFX miracevon; time for Raider or ASM to 
move this portion of calculated de- 
tection range. 


ISEED value. Initial seed value used in genera- 
tion of random numbers. Can be up to a seven 
digit number. 


X-array print out switch: 
0) =) nor print 
i=" print out X—-array as read in, target 
lists, and target list assignments. 


Random number counter. Set to O by model. 
Maintains running count of random numbers drawn. 
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X ( ) Inputs 


14 , Total number of random numbers used. Computed 
by model. 
15 Converted X-array print out switch: 


O = no print 
1 = print out X-array as converted by model. 


16 Blue first shoot control: 

O = Blue may fire as soon as Red detected. 

1 = Blue may only fire after first Red 
missile launch. 


17 Number of jamming raiders: 
no jammers 


1 = one jammer 
2 = two jammers, etc. 
18-19 Spares 
20 Altitude in feet above which a raider will be 


considered as being at "high altitude" for pur- 
poses of obtaining a visual identification of a 
non-radiating picket. 


ek. Probabulaty of “high altitude," from X(20), 
raider. launching ASM at a non-radiating picket. 


22 Probability of "low altitude," below X(20) in- 
put, raider launching ASM at non-radiating 
prcket. 

23 Maximum range (nautical miles) at which a raider 


Ts able to: obtain a2 visual identification of 2 
Blue surface unit. 


24 Geographical plots: 
0 = no plots requared 
1 = print plot of Blue positions 
2 = pring, plot Of rarder tracks 
3 = print’ both’ plots 
25-50 Spares 
(51-60) (Debug switches - cleared at end of first itera- 
ton aa pranvoute wal occur at onput value is 1.) 
51 Ship table. 
o2 Raid table. 
oS SAM table. 
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60 


(61-70) 


61 
62 
63 
64 
65 
66 
67 
68 
69 
70 
71-74 


75 


76 


77 


Inputs 


Random numbers used. 

Subroutine login/logout. 

SAM/ASM engagement, ASM impact, and SAM kill. 
ASM detection data. 

Closest point of approach and IBAT data. 
SAM-ASM Intercept data (angles). 

Store and Get Event Log. 


(Statistics generated by the model, ignore for 
input.) 


Total ASM's fired for all iterations. 
Total ASM's failed at launch. 

Total ASM's killed by SAM's. 

Total ASM's failed in flight. 

Total ASM's destroyed with launch platform. 
Total SAM launches. 

Total SAM failures. 

Total SAM overkills. 

Total ASM successfully launched. 
Total ASM hits. 

Spares 


Number of SAM rounds remaining in Fleet. Inter- 
nal use by the model. 


Hazard zone radius, nm. If predicted intercept 
OL tarcet aS within this ranee of another Ship, 
then the firing ship must hold fire until the 
predicted intercept point has moved beyond this 
range. 


Breakpoint altitude, ft, used to select proper 


SAM velocities, ranges, and intercept contours, 
tor “hieh os "low SAM contours. 
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mC) Inputs 





78 Ship position inputs coordinate system: 


O = Polar (000 True is Y axis). 
iv= Cartesian (xs=e ; Y= ). 
79-80 Spares 
81 Radarv horizon, factor): (default value 1.23). 
82 End of game time (default value 120 min). 
83 Launch: interval forsone raid, min. 
84 launch antervais fornaddi maids, min. 
85 Nunberr ot, shap corpwhiche the dasplay of hit fre- 


quency and hit interval is desired. Enter O if 
analysis is not desired or if the number of 
iterations is less than 10 or greater than 100, 
in which case the histogram can not be generated 
by the model. 

86 Ingexval for anatial Launch of Red SSM's, mani. 


87 Recochnition delay, in manutes, for stationary. 
Red raiders. For use when X(10) is not zero. 


88 Number of ships. 
89 Number of raiders (including Red aircraft which 
are only jammers, those without weapons, and Red 


surface/subsurface units).. 


90-100 Speres 


2. ohipeelass inputs 
The different types of BLUE ships in the simulation 
are specified by the ship class inputs. Included under one 
iiscing are all of the ships of a particular class (e.g. 
mG, CG, FFG) which possess Similar weapons and detection 
capabilities. A maximum of 10 different classes may be 
Mmeed-) Ship class 1 is input into X(101) to X(140); ship 


Brass 2 into X(141) to X(180); etc. 
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Te) Inputs 


101 SAM/SSM battery #1 class* 

102 Sam/SSM battery #2 class 

103 SAM/SSM battery #3 class 

104 SAM/SSM battery #4 class 

105 Search radar #1 class number” 

106 Search radar #2 class number 

107 Search radar #3 class number 

108 Search radar #1 antenna height (ft) 

109 Search radar #2 antenna height (ft) 

110 Search radar #3 antenna height (ft) 

et Fare control radar #1 antenna height (ft) 
a2 Fire control radar #2 antenna height (ft) 
113 Fire control radar #3 antenna height (ft) 
114 Fire control radar #4 antenna height (ft) 
15 Coordination capability for this ship class: 


O = none (override) 
1) = coordination, X(14), if» fleet. coordination 
TS. In, USE. 


116 Probability of correct identification and desig- 
nation to Weapon Direction Equipment (WDE). 


a7 Class of RTDF #2 (detection to ID and designa- 
tion to WDE). 


: Enter 1 if the battery has SAM/SSM's where class in- 
puts are listed in X(501) to X(540); enter 2 if SAM/SSM 
mipuLrs are in X(541) to X(580); ete. 


2 Enter 1 if the search radar class inputs for the re- 
Spective search radar are contained in X(901) to X(920); 
enver 2 if the: search radar: class; inputs are in X(921) to 
meo4+0); etc. © Note thateithe Fire Control Radar class is 
listed with its associated SAM class inputs. 
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124 


125 


126 


128 


129 


130 


132-140 


Inputs 


Class of RTDF #3 (designation to assign FCS). 
Class of RTDF #5 (firing delay). 


Class of RTDF #7 (decisional delay for horizon/ 
sector scan mode -- used with X(7)). 


Number of hits required to put ship out of action. 


Left SAM battery #1 launcher train angle limits 
(deg. true). 
Right SAM battery #1 launcher train angle limits 
(deg. true). 


Left SAM battery #2 launcher train angle limits 
(deg. true). 


Right SAM battery #2 launcher train angle limits 
(deg. true). 


Left SAM battery #3 launcher train angle limits 
(deg. true). 


Right SAM battery #3 launcher train angle limits 
(deg. true). 


Left SAM battery #4 launcher train angle limits 
(deg. true). 


Right SAM battery #4 launcher train angle limits 
(aee.™~ true’). 


Ability to determine target at which an ASM has 
been fired: 
ian0s= No 


Designation for non-radiating picket units using 
bistataceradar : 

O = No; unit will be targeted using only the 
normal, tarsetylist, algorithm. 

1 = Yes; unit will be targeted using both 
the normal target list algorithm and the 
additional algorithm for non-radiating 
pickets. 

If some bistatic units are to be used in the 
main force, a separate ship class must be en- 
tered for them. 


Spares 


50 







_ 








- 
= 


- 
_ 7 _ 
acti Jat inahy > OT 


tealeb caleft). oe age 
\AaSivial a4 <5 a \) Toe “Orie Ya = 


ony Tee, oe 
: nes OF 


i 
fl a 
boved 






os Ege 
° 


2 


X¢ ) Inputs 


141-180 ° Ship class 2 
181-220 Shapvelass’ 3 
221-260 Ship class 4 
261-300 Ship ‘class! 15 
301-340 Ship: class 6 
341-380 Shap class. 7 
381-420 Ship class” s 
421-460 Shap class 9 
461-500 pla Class 20 


3. SAM/SSM Class Inputs 


These inputs specify the BLUE SAM/SSM system capa- 
—pebities in terms of missile flight characteristics, fire 
control radar class, salvo firing doctrine, missiie relia- 
pwebitites and kill sneer ey A maximum of 10 different 


SAM/SSM classes may be input. SAM/SSM class 1 is input 


towo (501) to X(540); SAM/SSM class 2 into X(541) to X(580); 


etc. 

me?) Inputs 
501. Max range, nm, against "high" altitude target. 
502 Average speed, knots, for expected intercept 


range against "high" altitude target (converted 
by the model to nm/min). 





mC __/) Inputs 


530 ° A, nm, “high” altitude target” 

504 B, nm, "high" target” 

505 N; tom ’"hig¢h" target? 

506 SLOPE against "high" Creer 

507 Shik. hm, tor “haeh” target” 

508 Max range, nm against "low" altitude target 
509 Average speed, knots, against "low" altitude 


target (converted by the model to nm/min). 


BLO La nmagcarmst a “Low target 

opi ay nm, against a "low" target 

Ba NS for “tow "target 

513 SLOPE for "low" target 

514 SEEK, nm, for "low" target 

52:5 Minimum range, nm 

516 Maximum altitude, ft 

S17 Number of Fire Control Radar (FER) channels for 


FCR Type 1. (O=missile not requiring a dedi- 
cated fire control system. ) 


518 Class oF *FCR = Type 1 
519 Number of FCR channels for FCR Type 2 
520 Class of FCR - Type 2 
021 Number of FCR channels for FCR Type 3 
522 Class of ECR — Type 3 

3 


The definitions of A, B, N, SLOPE, and SEEK are con- 
tained in the description of subroutine TIMMY in Appendix B. 
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523 


537 


538-540 
541-580 


581-620 


861-900 


Inputs 
Horizon ysector,sean capability for fire control 
radars: 
O = No’ (used when X(7) set) 
1 = Yes 
input as XXX = unit digit FCR Type 1 
ten digit FCR Type 2 
hundred digit FCR Type 3 
Number of SAM's per salvo 
Launch cycle time, min 


Probability that launcher is available (used 
only when X(2) is set to 1). 


Magazine size 

Spare 

Probability of single-shot kill 

System reliability 

Reliability (launch, guidance, and fuse) 
Terminal guidance time, min (if HAW set to 0). 
Kill assessment time, min 

Spares 


Flag for Blue surface to surface missiles: 


Single shot kill probability for a crossing tar- 
get 


Spares 
SAM Class 2 


Sam Class 3 


SAM Class 10 
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4. Radar Class Inputs 
Search and fare ‘control radars for both BLUE and RED 


units are specified in this section of the inputs. Search 
radars classes are related to the ship class and raid class 
inputs. The fire control radars are used with the appro- 
priate SAM/SSM class inputs although their antenna heights 
ere contained in the ship class’ inputs. The radars are 
specified on the basis of their detection capabilities, beam 
widths, side-lobes, radar frequency, and delay functions. 

In the case of a bistatic radar, the location and antenna 
height of the transmitter must also be specified. A maximum 
of 15 radar classes may be included. Radar class 1 is input 


mmo ( 901) to XC920)> radar class 2 in X(921) to X(940); 


etc. 

x) Inputs 

901 Radar Type: 

O.= Buistatic 
1 = Monostatic 
eh 

902 B, detection range in clear environment, nm/m?* 
(see Appendix C). i 

903 a, detection,range with standoff jamming of 
lw/mhz, nm/m* (see Appendix C). 

904 Side-lobe ratio. Number of db that average 
Side-lobe is below main beam. Enter as a nega- 
tive number. Converted by model to side-lobe 
coefficient for use in obtaining jamming effects. 

905 Maximum instrumented range, nm. 

906 Class.of «reaction, timeydistribution function: 


Detectability-to-detection for search radars. 
FCS,assipned.,fo,acquisition for Lire,control 
radars. 
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mG 2) inputs 


907 Hadarasitrequency band dndexn.y Lgr2543,9Ore4. 
Used with jamming to provide for frequency dif- 
ferences in various radars. 


908 Cilass-of,reactaion time distribution: function for 
fire control radars: 
Detectability to detection delay in horizon/ 
sector scan mode. 
Blank for search radars. 


909 Max Director assignment range for FCR, nm. 
910 Spare 
(on i Radar beam width in degrees. (Converted to one- 


half beam width in radians by the model. ) 


oi Spare 

913 Probability radar is fully operable. (Used when 
KX(2) set.) 

914-916 Spares 

917 Ship number of bistatic transmitter. Used only 
Wit DiStatie radars. ff transmitter 1s. not 


located ona ship, "a “dummy “ship must be input. at 
the appropriate geographical location of the 
transmitter. 


918 Bistatic transmitter antenna height in feet. 
Used only with bistatic radars. 


919-920 Spares 


921-940 Radar “Class *2 


1181-1200 Hadar “Crass “l5 


So. lohip Position Inputs 
BLUE ship positions are input by specifying the 
Geass Of a Ship that is at a given location. The position 
may be entered as the range and bearing from force center or 


foe and Y coordinates’ in’ nautical miles. “The form must be 
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the same for all Bing a and is dependent on the value of 
x(7S8). The left and right firing bearing for the ships must 
aise be specified if firing sectors are used. If a mix of 
some ships with firing sectors and others without is used, 
X(8) should be set to 1 and appropriate bearings entered for 
all ships. Ships which are not limited should be assigned a 
fert..wearing of OOO and a right, bearing .of 359.99. A maxi-— 
mum otf 40 ships may be used. Ship 1 is input into X(1201) 


feerAsz05): Shap 2 in K(1206) to X(1210):. etc. 


X( ) Inputs 
1201 Ship class number. 
1202 LiIPKELS i= TOicaac lean nautieal males 
— i, & Coordinate in) nm 
1203 it- xC7S) = 0,, beariane from forces center in 
degrees true 
=, Y coordinate in nm 
1204 left sector Limit, degrees true; use if X(8) = 1. 
1205 Right sector limit, degrees true, use if X(8)= 1. 


In @ssignine sector limits to ships with no 
Limitations, put OOO mm X(1204) and 359.99 an 
x (2205:). 


1206-1210 Ship number 2 


1396-1400 Ship number 40 


6. Reaction Time Distribution Function Inputs 
these imputs provide appropriate delays, in minutes, 
for processing events during the simulation. A detailed 
description of the functions together with an example is 


contained in Appendix B under subroutine REACT. There are 


56 


ot as‘ > 


7 


Tea 2 
ee @ ae 

-_ . ; 7 an 
on eS £20 2:71 ifiea. bow Weroliees aed oe + 


% « 


bana of Kort iw aed > one ato ser Bee 





seven events which may be specified by up to five different 
classes of delay functions for each event. The reaction 
mameLions 2, 3, 5, and 7 are related to ship class inputs. 


Fumetions 1, 4, and 6 are used with radar class inputs. 


#1: Detactability to Detection (5 classes) 

a) Inputs 

1401 Delay in minutes for 0.1 cumulative probability. 
1402 Delay in minutes for 0.2 cumulative probability. 
1403 Delay in minutes for 0.3 cumulative probability. 
1404 Delay in minutes for 0.4 cumulative probability. 
1405 Delay in minutes for 0.5 cumulative probability. 
1406 Delay in minutes for 0.6 cumulative probability. 
1407 Delay in minutes for 0.7 cumulative probability. 
1468 Delay in minutes for 0.8 cumulative probability. 
1409 Delay in minutes for 0.9 cumulative probability. 
1410 Delay in minutes for 1.0 cumulative probability. 


1411-1420 Class 2 of-RIDFL 


1441-1450 Class Sof REDF1 


#2: Detection to ID and Designation to Weapons Direc- 
tion Equipment (5 classes) 


1451-1460 Class 1 of RTDF2 


1491-1500 Ciass 5 Of RIDE2 


57 





ec.) 
oN 


1501-1510 


(1541-1550 
#4: 


1551-1560 


1591-1600 
#5: 


1601-1610 


1651-1660 


1691-1700 


#7 


1701-1710 


1741-1750 


1751-1800 


inputs 


ID and Designation to Fire Control System (FCS) 


Assignment (5 classes) 


Class 1. of RIDES 


Class 5 of RIDES 
FCS Assigned to Acquisition (5 classes) 


Class i of REDE4 


Class 5 of RTDF4 
FCS Acquire to Fire (5 classes) 


Ciass ji tof VRLDES 


Class 5 of RIDF5 


Detectability to Detection for Fire Control Radar 


in Horizon/Sector Scan Mode (5 classes) 


Chass yi sof TRTDEG 


Class 3 Of RIDF6 


Decisional Delay for Horizon/Sector Scan Mode 


(5 classes) 


Class fof  -RepDr7 


Glass 65 cot IRL DER 


Spares 
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(6 pLSal Plses epmes 
Raid classes are used to ae the launch plat- 
forms for the RED ASM/SSM and standoff jammers. The launch 
Blattorms can be aircraft or surface units. Submarines are 
required to be on the surface for launch and are input as a 
surface unit. The raid classes contain the offensive capa- 
bility of the RED units in terms of weapons carried, radar 
class, jamming ability, and speed. The inputs also specify 
the radar cross section of the raider and the number of hits 
Torput it out of action. A maximum of 10 classes may be in- 
Meee hard class 1 is contained in X(1801) to X(1820); raid 


efass 2 in X(1821) to X(1840); etc. 


md) Inputs 
i80i Spare 
1802 Number of hits to put out of action, always 1 
for aircraft. Must be whole number. 
1803 Number of first type weapons carried. 
1804 Class of first type weapon. 
1805 Number of second type weapons carried. 
1806 Class of second type of weapon. 
1807 Radar class number. 
1808 Salvo switch (surface and subsurface only) 
O = 1 missile per salvo 
1 = 2 missiles per salvo 
1809 Launch, cyc¢le time for salvo, minutes for surface/ 
Subunits only. 
1810 Jamming power, w/MHz on band 1. 
1311 Jamming power, w/MHz on band 2. 
1812 Jamming power, w/MHz on band 38. 
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X( ) Inputs 
1813 Jamming power, w/MHz on band 4. 


1814 Raid speed, knots (converted by the model to 
nm/min). 

1615 Raid radar cross section, meters squared. 

1816 Fourth root of raid radar cross section (computed 


by the model). 


1817 Jammer beam width in degrees. (Converted to 
one-half beam width in radians by the model.) 


1818 Jammer side-lobe ratio. Number of db that aver- 
age Side-lobe is below main beam. Enter as 
negative number. (Converted by model to side- 
lobe coefficient for use in obtaining jamming 
effects. ) 

1819-1820 Spares 


1821-1840 Rard-Class 2 


1981-2000 Raita Class 10 


8. Raid Weapon Class Inputs 
The raid weapon classes specify the types of ASM/SSM 
by their flight characteristics, reliabilities, and radar 
eross section. A total of 20 different classes — be used. 
Weapon class 1 is input in X(2001) to X(2010); weapon class 


meen &( 2011) to X(2020); etc. 


X( ) Inputs 

2001 ASM/SSM speed, knots (converted by the model to 
nm/min). 

2002 ASM/SSM flight altitude, ft. 
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X( ) Inputs 


2003 Launch range’ 

2004 Reliabality (launch and initial guidance). 

2005 Inflight reliability (including homing and fus- 
ing). 

2006 Radar cross-section, meters squared. 


2007-2008 Spares 


2009 Weapon horizon, nm (computed by model). 
2010 Fourth root of radar cross Section Ceomputed by 
model). 


2011-2020 Weapon Class 2 


2192-2200 Weapon Class 20 


on Raid Position Inputs 
These inputs contain the location, course, and alti- 
puGe Of an individual launch platform. If the raider has a 
jammer, the jammer target is also specified. A maximum of 
S@ raiders is permitted. Raid 1 is input in X(2201) to 


meee); raid 2 in X(2207) to X(2212): etc. 


Bae.) Inputs 

2201 Rand class number (1,92, ete.) 

2202 Range, nm (from (0,0)). 

2203 Bearing, degrees true (from (0,0)). 


e Launch Range is measured from the position (0,0). The 
user must consider this fact when dealing with dispersed 
formations and multiple target list. The distance of the 
intended launch point from force center is the proper input, 
not the range of the launch point from the primary target. 
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X( ) Inputs 


2204 Altitude, feet (less than 50 feet for surface/ 
subsurface units). 


2205 Course, degrees true. 

2206 Jammer's jamming target: 
0 = no jam capability 
N = ship number of jamming target 
999 = force center (0,0) 


2207-2212 Raid number 2 


2734-2740 Raid number 90 


B. THE TARGET LIST ASSIGNMENTS 

The target list assignments are read in immediately fol- 
lowing the zero card indicating the end of the X vector 
inputs. Each RED raider must be assigned a target list, 
even if the raider carries no weapons. If the raider car- 
ries no weapons the target list information is not used, but 
the input logic requires that some number be entered. There 
must be a total of three cards which are to be dimensioned 
asot2. The first chery raiders are entered with card one, 
the second thirty raiders on two, etc. For example, if the 
Eamulation has four RED aircraft, column 2 card 1 would con- 
tain the number of the target list for raider number 1; 
eotumn 4 card 1 that for aircraft number 2; etc. Columns 9 
enrough 60 on the first card, and 1 through 60 on the next 


two cards could contain any integers desired, usually 00. 





©. THE TARGET’ LISTS 

The cards that input the target lists follow immediately 
after the three cards described above. The values to be 
read in are stored in the array PLIST. This #2 a (40,6) 
array in which the columns are the target lists. The array 
aetows Lor up to 40 individual targets and six separate 
histones of targets. The target list is identified by its 
column number, and it is this number that was entered in 
NaeGar.. The PLIST cards follow immediately: after the three 
input cards for NRTGT mentioned above. The format for the 
Petol cards is 10F8.4. The column Entries "are the total 
probability of targeting all ships up to and including the 


target ship in the given row. AS an example, in the case of 





Z 2 3 4 3) 6 
af 12 O20 29 0 O O 
2 42 29 500 O O O 
3 62 50 90 O O O 
4 82 eto £0 O O O 
5 1.0 0 £50 O O O 
6 O O 0 O 0 O 


The raiders that have been assigned target iis¢ tf would have 
feproObability of .12 of firing a missile at. ship 1 each time 
maares. it would have a probability of .3 of firing at 
moses, and .20 of firing at ship 3. <A raider firing a 
missile with target list 2 would have a probability of 0.0 
mraringe at ship 1 and .25 of firing at any of the other 


peur. A raider with list 3 would have.a .25 probability of 
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fecine at Ship 1, a .25 probability of firing at ship 2, and 
meso probability of firing at ship=4. 

The input cards for the PLIST array would be structured 
Poetolilows. Card one would contain .12-in columns 1-8, .42 
foe9-i16, .62 in 17-24, .82 in-25-32, 1.0 in 33-40, and zeros 
throughout the remainder. Cards two, three, and four would 
eenibain Zeros. Since card five starts list two, it is filled 
out in the same manner as card one, with, in this case, 
cards 6, 7, and 8 containing zeros. Since there are only 
mye Ships and three lists, the Last card required for this 
mun would be card 9. If the full 40 ships and six target 
lists were being used, a total of twenty-four cards would be 


mreequured) to input the target lists. 


64 


: : : . = Tt 


“4 7 sy = 


boetiecai te of Hlvoew vous Tera. wid 





. = yl 






PLD 


a 


7 ; 


i . - « t b 

cS. 6-2 enmwivoo. wi <!. aleites Sia eae 
*» 

—— ie YyrLkK OP=ER rk Pod : i--y Be Cm. iF fu 


"ee 
trys <4 , LS "eho 7 


VV... MODEL OUTPUT 


The model produces three types of output. The types are 
standard, special displays, and debugging. For debugging, 
there are a number of switches in the control section of the 
inputs which may be set to provide intermediary values dur- 
ing the iteration. The special displays provide one time 
outputs of geographical plots, histograms, and selected 
tables. The standard output occurs at the end of each itera- 
tion and contains an abbreviated battle history along with 


statistical data concerning the engagement. 


‘A. STANDARD OUTPUT 


ame standard output consists of the ASM table 


~ — Mee pey od 


, the SAM- 
ASM-RAID table, and statistics regarding hits, detections, 
and opportunities to engage. The contents of the tables are 
given below. The actual format of this output is provided 
in the tables of representative output at the end of this 
section. 
1. The ASM Table 
Column 1 ASM/SSM number. This number identifies 
the ASM/SSM by target number. The order 
in which the missiles were launched or 
destroyed can be determined by subtracting 
S90 from the target number. 
Column 2 Status code. This code indicates the 
disposition of’ the missile. The code is 
as follows: 


O'lDMil’ssile failed“on launch. 


Pismissrle was in “flaeht. 
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Column 


Column 


Column 
Column 
Column 


Column 


Column 


8 


9 


2 missile was destroyed by SAM. The 
SAM-ASM--RAID table will indicate which 
missile ship was responsible. 


missile failed in flight. 


co 


missile flight successful, missile 


impacted on target ship. 
5 missile destroyed when raider destroyed. 
Time of impact. This will be the time the 
missile impacted, or would have impacted, 


unless it was lost when raider was de- 
SEroyies, | in this case, 1 will be 0.0. 


Launch raid. 

Target ship. Target ship number unless 
the missile was lost with raider destruc- 
janie yay 

Launch time. The time of missile launch 
unless lost when raider destroyed. In 
thas case 2 1S the time of Voss. 

Weapon eeleassS.ey bE LOSt, with raid, 0... 


Mm coordinate at Jaunch. 


MeCoorginatre oat aunenh. 


The SAM-ASM-RAID Table 


Column 


Column 


Column 


Column 


Column 


al 


Target number. Targets 1 through 90 are 
raiders, 91 through 190 are missiles. 


Launcher. The launcher code AAB, where AA 
is the ship number and B is the launcher. 


Channel. The director-channel code AB 
where A is the director number and B is 
the channel number. 

Time Of. Launeh. Time of BLUE missile 
launch. 


Kill code. The code is as follows: a 
negative number, an overkill, the number 
indicating the number of missiles in the 
Salvo; a zero, a2 missile failure or system 
Support failure; a poSitive number, a 
kill, the valwe indicating the number of 
missiles in the salvo. 
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Column 6 Time channel free. The game time that 
this launcher-channel combination will be 
available for re-engagement or reassign- 
ment. 

oa, Ovabistical Output 
The output-form of the statistics for each iteration 
are underlined. The definitions of the terms used follow 
each form. 

The ships took "'H" hits: The total number of hits 

(H) sustained by the BLUE force. 

Ol ASiis tarred), UY" C12" pereent)) failed at launch. 

X: total number of ASM/SSM that survived to launch. 

Y:; of X, the number that failed on launch. : 

Di =. | Neopers. 

Of launched ASM, "W'"' ("P" percent) were killed by 

DaMs) 20) (Che Eper cen &) hailed hin Ghent: 

Launched ASM; from above, X-Y. 

W: the number of RED SSM/ASM shot down. 

P: W/Launched ASM. 

Qs lithe number (of RED (SSM/ASM that failed Gn flight. 

R: Q/Launched ASM. 

Of ‘total ASMs available "I" ("J" percent) were 

destroyed when mother killed by SAMs. 

Total ASMs: All SSM/ASM carried by engageable raid- 

ers or raiders that reach a launch 
position. 


I: The number of RED ASM/SSM lost when raider de- 
stroyed. 


J sh Totaly available, 
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"D" SAMs were launched of which "E" ("F'' percent ) 


failed) andy GiinG Wii mereent) were overkills. 


D: 


Ee: 


H: 


Total BLUE SSM/SAM launched. 


Number of SAM/SSM evaluated as failure at inter- 
cept eon sysrem SUPPOrt. 


E/D. 


Overkills. Overkill occurs when the target does 
not survive to the SAM intercept point due to 
beine Shot down. If the target is’ a missile and 
the RED missile reaches its target prior to 
intercept, this intercept is scored as a failure 
and not an overkill. 


G/ Di 


The ship statistics for each iteration are printed 


meatable labeled SHIP STATS. The format is as follows: 


Cotumn” 1 Ship number. 


Column 2 Chances to engage. The number of RED tar- 
gets that will enter this ship's missile 
envelope. 


Column 3 Number of engagements. The number of 


actual missile intercept evaluations. 


Column 4 Number of targets killed. 


Column 5 Bombers detectable. Of those-targets 


that will enter the ship's missile enve- 
lope, the number that are raiders. 


Column 6 Bombers detected. Of the raiders detect- 


able, those detected. 


Column 7 ASM's detectable. Those ASM/SSM that 


will enter the ship's missile envelope. 


Column 8 ASM's detected. 


The following additional statistics are provided 


wets the SHIP. STATS, table: 


Average detection range: The average detection 





range of each ship on those targets that 
entered its missile envelope. 
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Ux AsMsawere figed son ibystherfikkeet. xX is the num= 
ber of ASM/SSM that were engaged by one 
or more BLUE units. 


"X" Bombers were fired on by the fleet. xX is the 


number of RED raiders that were engaged 
by one or more BLUE units. 


ASMs fired at. The number of ASM/SSM that were first 
engaged by each ship. 


Bombers fired at. The number of RED raiders that 
were first engaged by each ship. 


A table of the frequency of hits received by ships 
as printed at the end of the run. This table displays a 
Srsteribution of the hit frequency over all ships for all 
fugs. Lt constitutes a histogram of hit frequencies for all 


Baum units. 


B. SPECIAL DISPLAY SWITCHES 

The below listed switches will provide the indicated 
SGutpurt if set to one or the specified value. Setting X(2) 
fo 1 will provide the equipment availability status for each 
iteration. The remaining switches provide a one-time output 
fom a run. 


X(2) Indicates equipment availability check desired. 
Prints out each ship number and equipment status. 


X(24) Plot of BLUE unit positions using DISPL if set 
toa. 


Prot pom raider, tracks and Jaunch points jusing 
RAIDIS if Set to 2. 


Both plots are provided if set to 3. 
mC51) Prints the SHIP table atimthe tend wf ithe first 
iteration. The table isastructured as follows: 


(see Table IV) 


Column 1 Ship number. 
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X(52 ) 


Column 2 Ship class number. 

Column 3 Ship's X coordinate. 

Column 4 Ship's Y coordinate. 

Column 5 Hits to sink.”’-The number of RED 
SOM/ASM missile hits to put out of 
actvon: 

Column 6 Class address. The ship class start- 
ing andex an) the XxX array. 

Column 7 Up/down. O = sunk 

1 = operating normally. 

Prints the RAID table®at~ the’ end of: the first 

iteration in the following form: (see Table IV) 

Column 1 RAIDER The raid number, i-90. 

Column’ 2 CLASS The class starting index in 
the X array. 

Column 3 TIMER (R) The time the raider was at 
the position indicated in 
columns 4 and 5. 

Column 4 XR The raider's X position at 
the time specified in 
column 3. 

Column 5 YR The sardéer’s ¥Y position at 
the time specified in 
column 3% 

Column 6 XEOG The X pesition of the 
raider at the time speci- 
fied as the end of the game. 

Column 7 YEOG The Y position at the end 
of the game. 

Column 8 XL The raider X position at 
its -tirst*Vauneh’ tame. If 
the raider has no missiles, 
or cannot reach a launch 
-position, the end of game 
position will also be indi- 
cated here. 

Column 9 Yi The raader Y position at 


rust awnen:: 
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X(53) 


Column 


Column 


Column 


Column 


Column 


Prints ythe 
iteration. 


Column 


O 


aATlumn 
CEM 


Column 


Column 


Column 


10 


al 


12 


13 


14 


a 


2 


ALT 


HORZ 


WEP-CLASS 


WEPS 


UP/DOWN 


The raider altitude in feet. 
If 150 vieethor "be Low jit ive 
raider is considered a sur- 
face t,or subsurface unit, 


The Prarder radar honrzon 
in miles. 


The class number of the 
raider's longest range 
ASM/SSM. 


The total number of weapons 
carried by the raider. 


The status of the raider 
at the end of the game. 
0 = ‘shot down or sunk 
1 = cperating normally. 


SAM table at the end of the first 
The table is organized as follows: 
(see Table V) 


SAM 

SAM CLASS 
INDEX 
SHIP 
NUMBER 


SAM LEFT 


CH1 FREE 


Wal 


The salvo number. 


The SAM tyve class number 
for the SAM/SSM in this 
BLUE miss-le salvo. 


The firi ne isin py 


The number of SAM/SSM re- 
maining in the battery's 

magazine after this salvo 
was fired. 


This is the time that chan- 
neti for the fire control 
radar used in this inter- 
cept will be free for re- 
engagement or reassignment. 
If the missile requires 
director illumination 
throughout missile flight, 
or the missile is totally 
independent of the direc- 
Lor Aethisecs Cant andi cavadon 
of system availability. 
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Columns 6-10 Channel availability times 
for systems employing a 
time sharing system of 
homing. 


Column 11 SALVO The number of missiles in 
the salvo. 


X(85) If the number of iterations is between 10 and 
100, prants, a hastogram of hit itrequency and 
ASM/SSM arrival interval for the BLUE ship who's 
number has been entered. 
C. THE DEBUGGING SWITCHES 

mae debugging switches, X-3,9,12:;19,.54,55,56,57,58,959, 
60, are set by assigning a value of 1 to the specified loca- 
tion in the X-array. These switches are subroutine related 
and will provide the indicated output whenever the program 
‘cycles through the appropriate portion of the subroutine. 
Reference should be made to the description of the subroutines 
in Appendix B when using the debugging switches. 

In the below format, the debug switch which must be set 
to 1 is listed and is followed by the subroutine name with 
ies associated output. The output is in the form as it ap- 
pears in the computer printout. The parenthesis indicate 
the location of the intermediary values. 

X-3 

SUBEL Detection Fange by ship ( ) om raider ( ) 
Wie cadarw() an clear, jamming, and 
horizon limited environments. 


Reaction time delays for raid ( *) by ship €°). 


Bustatic radar gain and lose contact times 
On raider © } by Ship: (_).. 
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X-12 


SUBE2 


PRIOR 


TIMMY 


MAIN 


ASM/SSM failed on launch. 


Detection range by ship ( ) on ASM/SSM ( ) 
With search radar ( ) in clear, jamming, 
and horizon limited environments. 


No detection due to reaction time delay. 


Detection range on ASM/SSM by ship ( ) on 
missile ( ) with system ( ) fire control 
nada (@ ). an “clear ,"jammine, and horizon 
limited cases. Jamming range of 999. indi- 
cates no jamming. 


Earliest search radar detection time and 
earliest fire control radar detection time 


CoC): 


When missile launch time will exceed the 
launch interval set by the user, the launch 
does not take place. Output is: bomber ( ) 
could not fire as launch time was ( ) and 
all ASM's had to be launched by ( ). 


EV 3 Ship ({-)--target- (€ ) reason ( ). 


EVs SAM (¢.) battery € ) channel € ). 

PRIOR debue data input channel ( ) output 
channel ( ) battery number ( ) input target 
@ ) eutput“tarcet (©). 

Intercept data: tarcet, firing ship, firing 
ship position, intercept position, number of 
firine Zones, lett and right bearings of the 
firing zones, launch bearing of' SAM. 


barce: farang ship, indicator flag, reason. 


Original X array as read into the computer. 
The target lists as read into the computer. 


The target list assignments as read into 
the computer. 
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EFFECT 
JULIE 
EVENT 
GETEV 
STOREV 
REMOVE 
TMDASM 
TMDBOM 
SUBE1 
SUBE2 
SUBE3 
SUBE4 
SUBE5 
SUBE6 
LOCTIM 
PRIOR 
TIMMY 
RCHECK 
JAMM 


BIDET 


The X-array as converted by the program. 


Prints out all random numbers used in 
Lies: Bheratnonm, 


used. 


Provides the computer time that the program 
enters and leaves the specified subroutine. 


Login, 
Login, 
Login, 
Login, 
Login, 
Login, 
Login, 
Login, 
Login; 
Login, 
Login, 
Login, 
Login, 
Login, 
Login, 
Login, 
Login, 
Login, 
Login, 


Login, 


logout 
logout 
logout 
logout 
logout 
logout 
logout 
logout 
logout 
logout 
logout 
logout 
logout 
logout 
logout 
logout 
logout 
logout 
logout 


logout 
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X-56 


SUBE4 EV 4 ASM ( ) ship ( ) ASM status ( ) hits 
ontishi ph (s )?tames (29% 


SUBES EV) Si target number’ (*) fringe ship'( ») 
Launcher ©); salvo. numbers (i )* tame (> ). 
X-57 
SUBE1 Detection ranges by ship ( )} on raider ( ) 


with system ( } fire controkt radar ( ) in 
clear, jamming, and horizon limited cases. 


Raider’ @) surface or submarine unit, ship 
( *{) has no SSM; 


Detectaon wanges' on raider ¢€ ) by ship ( ) 
in horizon search mode. 


SUBE2 Reaction time delays. 


Detection ranges in horizon search mode of 
ASM/SSM ( ) by ship ( ) in clear, jamming, 
and horizon limited cases. 


TMDASM No engagement 
skip ( ) batte 
ASM/SSM ( ) 


shows maid if \N toa 
- ta 
miss pile ship. -( 


No engagement possible; raider ( ) ship ( ). 


Engagement data —- time entering and time 
leaving a given range from ship ( ). 


X-58 
JULIE Number of SAM/SSM in each battery. 


SUBE2 Prints the NONO table after each raider's 
launch cycle. 


NONO table codes are contained in the de- 
seription of SUBEZ in Appendix 5. 


SUBES For saivo: ( ) results € ). 


us 
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X-59 


TIMMY 


X-60 
JULIE 
EVENT 


STOREV 


Crossing tarcet data: target, firing Ship, 
angle between target flight path and line 
of sight at launch, angle between SAM/SSM 
Pilaichtspath@and tarcet? flight path, angle 
between SAM/SSM flight path and line of 
Steht at Launch, target position ate launch, 
GSINE, range to target at launch. 


The events list at start of the iteration. 
The event retrieved. 
The event stored. 


The event lost. 
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VI. CONCLUSIONS 


The use of bistatic search radars is not a new concept. 
FADS provides a method of evaluating their potential useful- 
ness to the Navy. As demonstrated in the sample problem, 
FADS also has many applications in weapons systems analysis 
and the analysis of tactics. 

The model was structured sian the lines of the FLOATS 
model in use at the Johns Hopkins University Applied Physics 
Laboratory. The additional capabilities introduced by FADS 
should provide a basis for a more realistic view of the 
present day and future tactical environments. The employ- 
ment of decoys, dispersed formations, and passive missile 
traps are but a few of the ways in which FADS might be uti- 
jized. 

Lastiy, the model is ideal for use in an advanced war 
gaming course. The model is Bemisia NEE the simulation was 
written with ease of reading in mind. With rewriting the 
core requirements might be reduced, but, in its present form, 
FADS is easily read and understood, and it could be disected 


by the interested student with ease. 
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APPENDIX A 


BISTATIC RADAR SYSTEMS 


Although the bistatic radar system was described in the 
introduction as a radar system in which the transmitter and 
the receiver are not collocated it is necessary to point out 
that the separation between the transmitter and the receiver 
should. be comparable to the target distance, that is, on the 
order of several nautical miles as opposed to several feet. 
The assumed bistatic system is to provide detection and lo- 
e2aci10n information Similar to that obtained with current 
search radar systems. It should in essence provide the 
user with what would be on a plan position indicator (PPI) 
scope. No provision is assumed for height finding capabil- 
ity. The basic difference between the target geometry for 
monostatic and the bistatic radar systems is shown in figure 


ids. 





Figure Ala Figure Alb 


Monostatic Radar Principle Bistatic Radar Principle 





The pie eg ae iow: concept is not a recent development, 
although the technology required to bring the system up to 
the capabilities described above is. The bistatic concept 
predates the monostatic radar in that the development of the 
radar duplexer was required in order to convert the original 
bistaticyragar to-a monostatic.system. -Sholnik [Ref.. 10, 
pe 085] states: 


Its (the bistatic radar) principle was known 
and demonstrated many years before the development 
Omapractacal,monostatic.radarac,.2the,first.radar\ 
observations in both the United States and in 
Great Britain were made with Separated CW trans- 
mitters and receivers. These early radars were 
known as wave-interference equipment but were the 
same as what would now be called bistatic radar. 
Taylor and Young of the Naval Research Laboratory 
first demonstrated bistatic radar for the detection 
om shaps.i#n,1922 '—, Their, work. was,.disclosed, in, a 
patent issued in 1934. The early experiments with 
wave-interference (bistatic) radar led to the de- 
velopment of monostatic radar in the iate i930s in 
both this country and abroad. Further development 
was put aside after the demonstration of the more 
versatile monostatic-radar principle. Bistatic 
radar lay dormant for about fifteen years until it 
was "reinvented" in the early 1950s and received 
new interest. 


Current interest is conceivably motivated by considerations 
of electromagnetic counter-measures, defense against anti- 
radiation missiles, and deception roles. 

One method for target location with a bistatic system 
requires knowledge of the distance from the transmitter to 
the target, (RL), and to the receiver (RL); the horizontal 
angle between the base line (a line from the receiver to the 
transmitter) and the target bearing, and the elevation angle 
from the receiver to the target. However, in the simula- 


tion problem it is not necessary to determine a target's 
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Heocation since the target's position at any time is a func- 
tion of the input parameters. It is only necessary to de- 
termine if the target is detectable and, given that it is, 
determine where along the target's track an {aiyidual radar 
is able to detect the target. 

Detection occurs when the reflected power at the re- 
ceiver is greater than the detection threshold of the re- 
ceiver for a given probability of detection. The radar 
range equation for the bistatic radar is identical to the 
well known monostatic range equation except that the power 
received Var ie inversely with the product of BORD: Addi- 
tionally the two-way propagation losses differ since the 
-path of the reflected signal is not the same as that of the 
transmitted signal. The range equation can be readily 
written in terms of the constants for the radar system, the 
variables Ry» R.. and the target's radar cross section. 

For comparison purposes both the monostatic and the bi- 
Seat 1.c ferection range equations are given below: 


PG NGo 
. 4 a. hens tae 
Monostatic as (An) LER 


PuGaGs hac 


$ : = | Cael as 6 
SEOUEIS l= iain en henna 
Pp p s rx 


where: 
R = range to target (monostatic) 
R, = Trante from transmitters to Gareet C(bistatic ) 
R.. = range from receiver to target (bistatic) 
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Pp. = power of transmitter 

Pe " = power received at the receiver 

G = gain of monostatic antenna 

G, = gain of bistatic transmitter ant cana 

er = gain of bistatic receiver antenna 

r = radar wavelength 

Oo = target's radar cross section (see below for dis- 
cussion of bistatic and monostatic approximation) 

Ly, = one way propagation losses for monostatic 

L _(t) = propagation losses from transmitter to target 

Pp (Chastatic ) 

L _(r) = propagation losses from target to receiver (bi- 

Pp SEAtCLG) 

Le = system losses 


By defining Pa as the minimum power required by the re- 
ceiver in jecse i ete a given probability of detection 
with typical propagation and system losses and by grouping 
constants of the radar system, the above equations can be 
Sampiafied to read: ~- R= Bo? and ROR, = Bo where 8) is the 
range at which a particular radar system can detect a one 
square meter target in a non-jamming environment. Values of 
[eeor current U.S. Navy radars are listed in Ref. 9. 

When standoff electronic noise jamming is present, the 
range equations are defined in a similar manner: 


PiGo 
5 . 4 eat 
Monostatic: Rit 4mB,.(S/J)IP 


a als 
PGC? o 
we were 7 5 
Ne a gee ese 
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where: 


Reet burn through range monostatic 

B.. = receiver band width 

S/J = signal-to-jamming ratio (cmideseoeonacanaon 
losses ) 

YIP = sum of jamming power arriving at the receiver 


The jamming level from an individual jammer is obtained 
by multiplying the power of the jammer by the appropriate 
coefficient for the geometry of the radar and jammer main to 
Side-lobe ratios and dividing by the square of the distances 
between the jammer and the radar. The individual jamming 
levels are then added to obtain the sum of the jamming power 
geriving at the radar. 

As in the case with the non-jamming environment, the 
above equations can be simplified by defining a term a which 
contains the constants of the radar system and which is 
equal to the range at which a one square meter target can be 
detected in the presence of one watt of jamming power per 


megahertz. The monostatic range equation in jamming is: 


Pt 1 
= am (o/ HIE 


and the bistatic equation is 


ot, 
2 


ROR S-a(o/ZJP) 


In solving the bistatic detection problem in the simula- 
e1on algorithm, Ry and Ro have only the target's position in 
common. It is therefore necessary to state the product RR. 


as a function of the target's position and solve the 
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detection equation in terms of the time at which the target 
@an be detected. If a constant CHAR is defined to be the 
product of 8° times o for a clear environment and a? times 
o divided by the sum of the jamming power for the jamming 
situation, the detection equation then becomes ROR. = CHAR. 
The radar is able to detect a target whenever the range 
product squared is less than or equal to the constant CHAR. 
By expressing R.. and Ro. in .cerms of tame, the solution of 
the bistatic equation Re Ge BCE) = CHAR = 0. provides the 
seames at which detection can occur and then be lost as the 
target moves along its track relative to the transmitter and 
the receiver. The below figure is a layout of the track and 
detection points that is illustrative of the detection equa- 


tion: 





Figure A2 


Bistatic Detection Geometry 
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Labels are as follows: 


Receiver 
XR 


YR 


Transmitter 


XX 


YX 


ZXN 


XXC 


YXC 


CPX 


Value 

& ‘Coord of posi tivon 

Y Coord, of posit ron 

antenna height in nautical miles 
x ‘Coord *of track CPA 

Y *Coord “of “track "CPA 


Distance to CPA 


ini Fea l posi eLon 

altitude in nautical miles 
terminal position 

velocity in nm/min 

velocity -in Xedireétion 
Velocity ian Y°direction 

X Coord of detection position 
= (X1+VxT) where T is the time 


of detection 


Y Coord of detection position 
= (YisVyT) ; 


The equation Ry(t)RE(t) = CHAR. =.0 is expressed in terms 


of the stated positions; time is the independent variable. 


For example: 


Re (¢) 


(XR-XRC )*+(YR-YRC )?+(ZRN-Z1)?+(XRC-(X1+VxT) )? 


+(YRC-(Y1+VyT))? 


(Vx? +Vy” )T?-2(Vx(XRC-X1)+Vy (YRC-Y1) )T+(CPR)? 


eCRRC-X1y PC YRC-¥1 9 


AT?+DT+E 


&9 





In a Similar manner RY (t) = AT*+BT+C where B = -2(Vx(XXC-X1) 


Avy YXC-Y1))T-and C = (OPK )C + CXC- 1) = CY XC-VN)* 
mae range equation in terms of t then reduces: to 


BELCD 4. CH = CHAR 
Se ey 


r+(2Pyr3+(S* 


A 





+ ao yt? +( 


ihe real roots of this Tourth-order polynomial in ascend— 
ae, order are the times at which the target is first detect- 
able, lost, redetected and then lost again. In the case 
where there are no real roots, the target is not detectabie 
Monts track. lf-only one-pair of real roots exists, the 
radar has only one detection zone against the target. 

The multiple detection zones occur when the target pene- 


ale 
trates the contour of the equation R i = CHAR? in more than 


™ 


two Jocations. These contours plot as ovals of Cassini in 
the plane containing the base line between the transmitter 
mie Lece1 ver... The contours j (shown) im figure As, repre— 
Sent curves of constant sensitivity for the bistatic radar 
system and a given target. In the case of low sensitivity, 
the curves consist of two disjunct ellipses with the trans- 
mitter and receiver at the furthest foci. As sensitivity is 
increased the curves become ovals about the receiver and the 
transmitter and eventually become an ellipse with the trans- 
mitter and receiver at the foci. If this ellipse is rotated 


about its major axis, a prolate spheroid is generated which 


Gefines the detection surface for the bistatic radar. 


90 








Figure A3 


Surves of Constant Sensitivity for Bistatic Radar System 


In addition to the requirement that a solution to the 
radar range equation exists, the target must also be above 
both the transmitter and the receiver radar horizon for a 
Betection to occur. The radar horizon restrictions can be 
seen in figure A4, which uses the same track as presented in 
figure A2, along with an additional axis representing the 
asiceratt position with respect to time. 

The positions in figure A4 are defined in figure A2. 
HORX and HORR are the distances to the transmitter and re- 
ceiver radar horizons for their respective antenna heights 
ema the target's altitude. TR1 is the time the target first 
enters the radar horizon of the Feceivier, Ikh2 tS the’ time £L 
aeparts. “TX and TX2-are Similar’ times for the transmitter. 
The time that the target is simultaneously above both hori- 
zons is the interval Ae where ci = max (TR1,TX1) and 
i = min (TR2,TX2). 
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HORX 


Figure A4 


Intersection of the Bistatic Radar Horizons 


The intersection of the joint horizon time interval 
fet) with the time intervals obtained from the range de- 
tection equation provides the time of detection and the time 
Or Jost detection for a target. 

These times, in conjunction with the ranges which are 
related to the times, are then used in the simulation to de- 
termine “shoot events" for the defensive forces and to accu- 


mulate statistics on the detection ranges. 
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APPENDIX B 


DESCRIPTION OF THE PROGRAM 


This appendix provides a description of the FADS com- 


puter program. The MAIN portion of the program is presented 


mars. . 


The principal subroutines follow with the minor sub- 


mMOouULInes and functions covered last. 


te OVERVIEW 


This paragraph gives an overview of the functions of 


the individual sections of the program. 


MAIN 


EFFECT 


JULIE 


EVENT 


REMOVE 


TARGET 


SUBE1L 


SUBE2 


SUBES 


The MAIN program reads in the input values, com- 
piles statistics and controls the simulation through 
the specified number of iterations. 


Subroutine EFFECT controls the individual itera- 
tion by selecting the proper subroutine to play 
the next event in the events list. 


subroutine JULIE initializes or resets values for 
each iteration. 


Subroutine EVENT controls the input and retrieval 
of events from the events list. 


Subroutine REMOVE retrieves higher priority tar- 
gets from the events list. 


Subroutine TARGET returns all of the "shoot" events 
which occur at a specified time. 


Subroutine SUBE1 computes detections of the raiders 
by the BLUE forces and determines ''shoot" events 
against the raiders. 


Subroutine SUBE2 computes the RED ASM/SSM launches 
and subsequent detections followed by the "shoot" 
events against the ASM/SSMs. 


Subroutine SUBE3 evaluates the engageability of a 


target based on the situation at the proposed BLUE 
missile launch time. 
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SUBE4 


SUBES 


SUBE6 


TIMMY 


JAMM 


BIDET 


DTIMD 


CPA 


REACT 


PRIOR 


RCHECK 


TMDASM 


TMDBOM 


ZLNKA 


TJSORT 


Subroutine SUBl4 determines ASM/SSM impacts on 
BLUE ships and evaluates out of action status for 
them. 


Subroutine SUBES evaluates the BLUE missile inter- 
cepts. 


Subroutine SUBE6 releases the fire control radar 
at the time of intercept. 


Subroutine TIMMY computes the predicted intercept 
point and flight time for the BLUE SAM/SSM and 
determines whether the intercept is within the 
missvle*s°capabwlaty: 


Subroutine JAMM determines the range at which a 
monostatic radar can detect a target in the pres- 
ence of jamming. 


Subroutine BIDET computes detection ranges and 
times of detection for bistatic radars in a clear 
or jamming environment. 


Subroutine DTIMD determines the real and imaginary 
roots of the fourth order polynomial used to solve 
for the bistatic radar detection times. 


Subroutine CPA computes the coordinates of the 
efosest “point of approach Vof “a’-target ‘to’a-unit. 


Subroutine REACT computes reaction time delays. 


Subroutine PRIOR determines the most threatening 
target among those that a defensive system may en- 
gage at a given time. 


Subroutine RCHECK evaluates a prospective target 
ship to determine if the target is above the 
ravder *s “radar horizon and within the “raiders 
maximum theoretical radar range. 


Subroutine TMDASM determines if and when a RED 
target will close within a specified range of a 
BLUE unit. 

Subroutine TMDBOM determines the time when a RED 
raider will come within a given range of force 
center. 

Subroutine ZLNKA is a sorting routine. 


subroutine TISORT as a sorting routine. 
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LOCTIM - Subroutine LOCTIM computes the raider position at 
a given time. 


ASMTIM - Subroutine ASMTIM computes the ASM/SSM position at 
a given time. 

RANDOM -—- Subroutine RANDOM generates uniform (0,1) random 
numbers. 

HISTG - Subroutine HISTG provides a histogram and asso- 


eaated stati sticsootshitspuhitanintervals) \andahats 
per minute for a specified BLUE unit. 


DISPL - Subroutine DISPL provides a geographical plot of 
the BLUE unit positions. 


RADIAN - Function RADIAN converts angles in degrees to 
angles in radians. 


RAIDIS - Subroutine RAIDIS provides a geographical plot of 
the RED rarder tracks and launch points. 


SETIMX - Subroutine SETIMX provides the internal computer 
clock time. 


ho 


MAIN 
The simulation is initialized by MAIN which also com- 
piles iteration results for the summaries. The MAIN program 
begins the first of its three principle functions by zeroing 
specified setae, It then proceeds to read in the X array, 
the raider target list assignments (NRTGT), and the target 
fists (PLIST). The second principle function is to enter 
moose initial values that will remain constant through all 
meeracions into the various tables. MAIN's third principle 
funetion is the accumulation and printing of the statistical 
results. 

The initial step in performing the above functions is 
the zeroing of the X array. In the next step, data are read 


mite the Xrarrayliby reading data cards in the format 
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© ,5F10.3)... The aries Vein is first read into an integer, 
ment, and the array RXATA(5). From ubeue. using IRT as the 
initial X address, MAIN fills the next five X array slots 
sequentially. xX array addresses that are to remain 0.0 need 
not be read in unless they are within the neighborhood of 
the five addresses contained on an input card. In this case 
they must be specified 0.0. The sequence of reading data 
Sards and filling the X array will continue until a card is 
read with the value O for IRT and 0.0 for the remainder of 
moe. card input. This is the flag for MAIN to proceed to the 
next step. 

The program then obtains from the values previously 
read into the X array the starting values for ISEED, the 
seed value for the random number generator, and NI, the num- 
ber of iterations to be played. End of game time is deter- 
mined and the X array is printed as it was read in. MAIN 


now computes the weapon (RED ASM/SSM) horizon as follows: 
a 
HORZ = (CURVATURE FACTOR) x (CRUISING ALT)? 


HORZ iS in miles and ALT is in feet. MAIN next converts all 
courses and bearings from degrees true to degrees in the 


cartesian coordinate system using the relationship: 
DEG(C) = 450.-DEG(T) 


Ali velocities are converted from knots to nautical miles 
per minute. Jammer and radar beam widths are converted to 
half beam widths in radians, and MAIN finds the fourth root 


of the raider and RED weapon radar cross sections. Next the 


96 


el 
iy - Be Ai Bank “Cass eee 
MNRIN I He DIDI Spon, tee 


a 
e ~-10 Pa ; teh <9 oy - 


wit ae TAY Sufeu, ,esils mp 
: id a” le 4 


‘ 
i 





target list assignments and target lists are read in, and, 
if X(12)is set, MAIN prints the X-array, PLIST, and NTRGT as 
read in. 

MAIN next fills out the ship table and the raid table. 
It counts up the number of BLUE missiles (both SAM and SSM) 
ayeatable, and, if X(15) has been set, MAIN prints the con- 
wemted X array. If X(55) is set, MAIN calls SETIMX to start 
meee internal timer. 

MAUN now calls EFFECT to control the first iteration. 
Upon completion of the first iteration, MAIN prints the ASM 
table and the SAM-ASM-RAIDER table. It then computes and 
Drants various statistics for this iteration and computes 
fad eprints cumulative Statistics for ‘all iterations thus far 
eemplteted. If X(51) is.set, the SHIP table is printed at 
the end of iteration one. If X(52) is set the RAID table is 
Printed and if X(53) is set the SAM table is printed. If 
X(54) is set, all the random numbers used in the first iter- 
auton are printed, in the order of use. 

Bit, the end iof iteration’ one," MAIN’ sets flags 9, 52, 53, 
memeso. Of, OS, 59, and 6GO to 0. MAIN then loops back and 
calls EFFECT to begin the next iteration. Upon completion 
ePethe last iteration, additional statistics are printed 


and the program stops. 


=a SUBROUTINE EFFECT 
HPERCT controls cach iteration. If X(55) is set for 


debugging, EFFECT calls subroutine GETIMX on entry and exit 
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to login and logout. EFFECT next calls JULIE to set the 
tables and arrays that were altered during the previous 
meration back to their, proper starting values. The simula- 
tion begins when EFFECT calls subroutine GETEV for the first 
event to be played. GETEV supplies EFFECT with the type of 
event to be played, NEXTEV, the time of the event, and a 
code identifying the type of event. If the event type re- 
turned by GETEV is O, the iteration is over and the sub- 
maine exits. If tthe event type is greater than, 6, the 
program prints an error message and terminates. Otherwise 


EFFECT calls one of the following subroutines to continue 


the play: 
SUBE1- for a RED raider detection and intercept 
ev elit,. 
SUBE2 for a RED missile launch-target-—detect- 


intercept event. 
SUBE3 for a BLUE missile firing event. 
SUBE4 for a RED ASM/SSM impact event. 
SUBES5 for a BLUE missile kill event. 


SUBEG6 froma “set fLirescontrods channel free" 
event. 


Upon return from one of the SUBE's, EFFECT again calls GETEV 
for the next event to be played. Event types used elsewhere 
in this manuscript correspond to the SUBE subroutine which 


plays the event. 


4. SUBROUTINE JULIE 
Subreutine JULIE initializes: or.resets those items that 


are required for each iteration. ft, C55) asset, JULIE 
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calls. GETIM upon entry and exit. JULIE is called by EFFECT 
and calls RANDOM, TMDBOM, STOREV, and EVENT. 

EVENT is called to zero the events idst;, andstorenter 
the end of game event into the events list. The following 
counters are then set to zero; NASM, NSAM, IFAG, IFBG, ITEMS, 
IFLASK, and IFIRST. The following arrays are zeroed; JAM, 
MiaT, ISAG, ISBG, STATS, TGTAV, IGAGE, NONO, JACK, and 
SAMASM. The RAID table and SHIP table are reset and the SAM 
table is zeroed. The missile batteries are reloaded and, if 
X(2) is set, the availability of each BLUE radar and launcher 
is determined. JULIE does this as follows. The first check 
is on the search radars; if at least one is operating, JULIE 
determines which fire control radars are operable. If no 
Search radar is available, the program checks tc see whether 
the horizon search mode for fire control radars is author- 
ized, X(7), and if so, whether the ship in question has 
radars with this capability. If the ship has a search radar 
and fire control radar combination, or a search radar com- 
bined with active, HOJ, ARM, or some other type missile 
which does not require a fire control radar, the launch 
availability is checked. If the search radars are inopera- 
tive, the fire control radars are horizon search capable, 
and at least one is operable, the launchers are checked for 
availability. If the ship has no means of firing or of de- 
eect ron) “it is marked 2 non-combatant anduthe next unit is 
evaluated. If the launcher can be coupled with a search/ 


detection/acquisition system, the ship is marked as a G-ship, 


79 


a ars . 
at be - 7 


1) a - 7 - ' > : 
TRAST yd boklesser She See 
ie 2 


; TVG Lee CE ae 


= 


s T5W ii hi a] vp’ 





the equipment status is entered in the NADAR array, and the 
next unit is evaluated. 

For each RED raider, JULIE calls STOREV to store a 
type 1 event at time O., and, if the raider carries ASMs or 
SSMs, a type 2 event is stored for launch time as determined 
by TMDBOM. 

Finally, JULIE totals the number of BLUE SAM(SSM) bat- 
teries available, sets the RELOAD array to zero, and, if 


mCGO) “1s set, prints the events list. 


5. SUBROUTINE EVENT 

Subroutine EVENT is called by JULIE. It enters the 
end of game time in the events list. If X(55) is set, 
EVENT calls GETIMX upon entry and exit. The event counter 
LOVER: is set to one and the events list. is zeroed. End of 


game time is set equal to the integer IT using the equation: 
Eta ehO00 26 oC Toads O-. 0005) 


muerte IT is equal to the value stored in X(82). ‘The end of 
eeme event, a type © event, is then entered in the events 
fist along with the time IT. EVENT then returns to JULIE. 
a. Entry GETEV 

GETEV is an entry into Subroutine EVENT. It is 
eatled by EFFECT and calls GETIMX at entry and exit if X(55) 
is set. GETEV removes the first entry in the events array, 
MNVENTS (1,1) and IVENTS (1,2), ‘converts these back into 
event time, type, and participating units, and sends these 


three. pieces of information back to EFFECT. The counter 
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LOVER is decremented and the entire events list is moved up 
one position. If the debugging switch X(60) is set, GETEV 
prints out TIME, NEXTEV, and IUNIT. 

ba > Entry. STOREV. 

STOREV is an entry in Subroutine EVENT and its pur- 
mose is to place the given event in its proper place in the 
events list. STOREV is called by JULIE, SUBE1, SUBE2, SUBE3, 
BUBES, and PRIOR. If X(55) is set, STOREV calls GETIMX upon 
enury and exit and prints the login/logout time. The inputs 
to STOREV are the time of the event, the type of event, and 
eae code for the units involved. STOREV first checks to see 
whether the time the event is to be played exceeds the 
allotted game time. If it does, STOREV prints a message 
fiving the time of the last event, tyne of event, and plaver 
code. | If LOVER is greater than or equal to 1000, the event 
is lost and STOREV prints a similar message. If both tests 
are passed, STOREV combines the event type and player code 
antLO one 9 digit integer (ICODE). Time is converted to an 
integer using the same equation as is in EVENT. STOREV 
places the event in the events array by time. If the time 
is not unique, it proceeds to compare the stored values of 
ICODE to the one to be entered. These are ordered from 
smallest to largest. When the proper place is found, LOVER 
is incremented, and the: remainder of the array is moved up 
ene address. If no place is found to insert the event, in- 
cluding even just right before the end of game event, an 
error message is printed and the game stops. If the debug 


ewatch X(60) 1S set, STOREV prints the event stored. 
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6. SUBROUTINE REMOVE 

Baum Tae REMOVE is called by PRIOR when PRIOR dis- 
covers a more threatening target for the missile system 
being evaluated. If X(55) is set, REMOVE calls GETIMX upon 
entry and exit. REMOVE receives the code IUNYT from PRIOR. 
IUNYT is the (event type/player) combination that is stored 
in the events array. IUNYT was found by TARGET and evalu- 
ated by PRIOR as the next event to be played. REMOVE 
searches the events array for a matching code group. If it 
finds the group it removes both the code and the time and 
then moves the remainder of the array up one, and exits. If 
REMOVE can not find the code, it prints an error message. 


REMOVE then prints the events array and exits. 


a. SUBROUTINE TARGET 

Subroutine TARGET is called by PRIOR. Its purpose is 
to provide PRIOR with all Type 3 events which occur at the 
present game time. The input for TARGET is time and the 
outputs are IUNIT and? MORE. “St clr IRS Gs ‘equal tovzero, 
TARGET sets it equal to one. TARGET then converts time to 
the integer format of EVENT and compares it to the time in 
Sie array location indicated by IFIRST. If.the times are 
equal, TARGET checks the event code. If it is a Type 3 
event, TARGET increments IFIRST, sets MORE equal to one, and 
exits with the event eodee If the times are not equal, or 
the code is not Type 3, MORE and IFIRST are set to zero and 


the subroutine exits. If MORE was set to one, PRIOR will 
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penivagain until MORE is set to zero» If X(55)°is set, 


TARGET calls GETIMX each time it is entered or is exited. 


8. SUBROUTINE SUBE1 

SUBE1 is the RAIDER detection subroutine. It is called 
by EFFECT and in turn calls subroutines CPA, JAMM, STOREV, 
REACT, EOCTIM,- RANDOM, IMDASM, BIDET, and, if X€55) is set, 
GETIMX on entry and exit. SUBE1 tracks the raider from its 
start of the game position to its end of game position and 
determines which BLUE units, if any, will be able to engage 
the raider. If it is determined an engagement is possible, 
detection delays, command and control delays, and equipment 
availability are checked to determine whether or not to 
enter a Type 3, or shoot, event in the events list array. 

SUBE1 first sets flags NASA, NSURF, and NJ. If X(16) 
is set, SUBE1 checks the ASM counter to see if the BLUE 
force has been fired upon. If X(16) is set and NASM is zero, 
SUBE1 stores a Type 1 event at game time plus 2 minutes and 
then exits. Otherwise, the raider's altitude is checked to 
determine whether the raider is an aircraft or surface/sub- 
Surface unit, and then the applicable raider characteristics 
are drawn from the X-array. 

Kach BLUE ship is next queried to determine whether it 
Bem engage the raider. If the ship has been sunk, or if it 
is a non-missile ship, the NONO table is set to 7 and the 
next unit is evaluated. The maximum missile range of the 
Operable systems on the BLUE unit is obtained. The dis- 


tinction between surface and air targets is made by comparing 
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target altitude to the fifty-foot ceiling for enemy surface 
units and the fifty-foot floor for enemy aircraft. A mis- 
sile system must be designated as a surface to surface, or 
surface to air system when the system characteristics are 
Zputc. there. is no cross matching of systems to targets. 
SAM's are fired only at aircraft and ASM's, SSM's only at 
ships. 

TMDASM is called to determine whether the target will 
enter the BLUE unit's missile envelope. If not, the NONO 
table is set to 9 and the next unit evaluated. If the tar- 
get will enter the envelope after the end of game time, the 
NONO table is set to 8. 

When it has been determined that the target can be en- 
gaged, the target CPA is obtained and the ship's engagement 
opportunity counter is incremented. 

The detection phase begins with a check of the search 
radars. Each BLUE unit is allowed up to three. If the 
radar is operable, the clear detection range, jamming envi- 
ronment detection range, combined radar horizon, and theo- 
retical maximum radar detection range are compared. The 
minimum of these is the radar's detection range for the tar- 
Bet an question.| This range is compared to that of the other 
search radars on board, and the maximum is taken as the 
Ship's detection range. 

The between radar comparison is made by converting de- 
tection ranges into time and adding in the reaction time 


delays obtained from REACT. 
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ifs there) is - search radar operational, SUBE1 checks 
to see whether the target can be etwewed in the horizon 
search mode. If there has been no detection and the BLUE 
unit is operating a bistatic search radar, the subroutine 
checks to see whether BIDET indicated two detection zones. 
If so, the search radar portion of SUBE1 is redone using the 
second bistatic zone. If there remains no detection, or no 
radar capable of detection, the NONO table is set to 2 and 
the next BLUE unit evaluated. 

Given a search radar detection, SUBE1 proceeds to eval- 
uate the missile systems for fire control acquisition. This 
is done by obtaining the maximum fire control radar range 
from among those radars available, and converting this toa 
time as was done in the search radar algorithm, If the mis- 
Sile of a system requires no fire control radar, the maximum 
missile range is used as that system's acquisition range. 

Having determined the maximum acquisition range, the 
probability of correct identification is compared to a ran- 
dom number drawn from RANDOM. If the ID check declares the 
target non-hostile, the NONO table is set equal to 3 and the 
next BLUE unit is evaluated. Otherwise, the time the target 
will be acquired is determined by calling TMDASM. 

The reaction times for weapons direction equipment de- 
lays, target assignment, and acquisition (if the missile re- 
quires a fire control radar) are obtained from REACT and 
added to the search radar detection time. The last time to 
Shoot is determined using the time the target leaves the 


maximum fire control range, and the SAM/SSM velocity. 
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The firing delay is obtained from REACT and added to 
the eater of the detection time and the acquisition time. 
The "shoot" time is compared to the” last time to fire, and 
af the shot is still possible, the program stores a shoot, 
Type 3, event for each director-launcher combination in op- 
eration. If the missile requires no director, the number 7 
is stored as the channel number, and only one shoot event is 
stored for that system. If the unit requires a director, no 
channel number is assigned. 

If the delays have precluded a firing, SUBE1 checks to 
see whether there was a second bistatic detection. If so, 
it returns to the search radar portion of the algorithm and 
proceeds as above evaluating the second bistatic detection 
zone. Otherwise it sets the NONO table equal to 5 and ynroc- 
esses the next BLUE unit. 

The horizon search portion of SUBE1 is entered as de- 
scribed above. The program checks the NADAR table to see if 
the directors can operate in the horizon search mode. If 
so, the program proceeds in the same way as the search al- 
gorithm, using the fire control radar parameters. Having 
determined the detection time, the decisional and reaction 
delays for horizon scan mode are obtained from REACT. The 
shoot time is determined to be the detection time plus the 
acquisition time, plus the maximum of the firing time and 
décision to fire time. A range check is made by calling 
TMDASM and if the time to fire is less than the last time to 
fire, a Type 3 event is stored for each director operable in 


the horizon search mode. 
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Atter all the BLUE units have been processed, SUBE1 


exutcs calling, GETIMX if flag X(55) has been set. 


2. SUBROUTINE SUBE2 

SUBE2 is the ASM/SSM launch and detection subroutine. 
fis Called by EFFECT and in turn calls RCHECK, REACT, 
STOREV, BIDET, RANDOM, TMDASM, TJSORT, CPA, LOCTIM, JAMM, 
mma at X(55) is set, GETIMX on entry and exit. 

The inputs to SUBE2 are launch time and raider number. 
The subroutine first checks the RAID table to see whether 
the raider Hale already been shot down or sunk. If so, the 
subroutine exits. Next the program checks the flag IFLASK 
to determine whether any RED ASM/SSM's have been fired. If 
IFLASK is greater than zero, the subroutine proceeds to 
check whether the launch time is greater than TMAX and exits 
meat is. If IFLASK equals zero, SUBE2 sets IFLASK to one 
and TMAX, the time by which all ASM's must be launched, to 
present time plus X(84). SUBE2 now retrieves the raider 
speed and, if the raider is an aircraft, determines the ASM 
launch interval, sets IFLASK equal to 2, and determines the 
total number of ASM's on board. If the raider is a surface 
unit; and no ASM's have been fired, IFLASK is set to 0. 

SUBE2 determines which of the two possible missile types 
has the longest firing range, (firing range is measured from 
coordinate (0,0)). The missile parameters are drawn from the 
X vector, and the launching sequence begins. 

The launching of ASM's is done by drawing a series of 


uniform (0,1) random numbers, depending upon how many missiles 
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are to be launched, and ordering them. .The exact launch 
time of each ASM becomes the present time, plus the next 
random number times the launching interval. The launching 
of SSM's from surface units begins with the finee launch at 
time plus some portion of the interval specified in X(86). © 
Subsequent surface launches are at fixed increments, X(1809). 
Having determined when, and in the case of surface to 
surface launches, how many RED missiles are to be launched, 
the targeting phase beings. LOCTIM is called to insert the 
maaaer poSition at launch into the RAID table. The raider 
target list is located and a temporary target list is con- 
Structed. Each ship on the raider's target list is cycled 
. through RCHECK to determine whether it is afloat, above the 
radar horizon and within radar faneee If all these tests 
are passed, the ship is entered on the temporary target list 
Byepucting the values from the original-target list in the 
next slot of the temporary target list. Once the target 
fist’ has Bedi formed, a random number, Z, is drawn from 
RANDOM and the first ship with a cumulative probability of 
being targeted greater than Z is designated the target. If 
there is none, the first ship on the list is taken. A check 
is made to see whether this ship has been designated a non- 
radiating picket. If not, SUBE2 proceeds to the detection 
phase. If the target chosen is a non-radiating picket, 
checks are made to see whether the raider can correctly 
identify the target. This is doen by determining whether 


ene raider: is within’ X(23) miles of the target. This 
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measurement is done in the horizontal plane. If the raider 
is close ‘enough to the target, RANDOM is called for a random 
number to compare with the probability of correct identifi- 
Car2on. lf the picket fails to qualify as a target, Lt is 
eliminated from the target list and the process starts again. 
Tf, at any time the situation arises where there are no tar- 
Bere On a raider's target list, launching stops and an event 
Type 2 is stored for time plus three minutes. 

Upon acquiring a target, the ASM/SSM is considered 
launched. A random number is drawn and compared to the 
fauneh reliability. The ASM table is filled out and, if 
launch was unsuccessful, the NONO table is set equal to 1 
and the next missile processed. Upon an unsuccessful launch, 
a check is made to see whether there is to be a second SSM 
in the same salvo. If so, the time is not advanced for 
this launch. 

Assuming a successful launch, the detection phase be- 
gins. The target ship's position is stored in the RAID 
table for use in BIDET, but otherwise, the detection se- 
quence follows that outlined in SUBE1. 

When all of the first type of missiles have been launched, 
and if the raider is carrying a second type of missile, 

SUBE2 decides whether or not the raider has reached its sec- 
ond launch point. If so, the appropriate parameters are ex- 
tracted and the launching continues; otherwise an event Type 
Eats Stored for the new launch time and the subroutine exits. 
When there are no second type missiles, the subroutine exits 


aitter the last launch. 


Oo 







oot. 7 
aT _ 


Wahine. b Wor hufkaotet 


aS lwiht Yostaron bo ext 




















- v 
on af si 


iS 


) uit Oa Sell Sosa i 





i] af i 
if ms sat 5TH ive 
iMrTyiasy 7” 
: . 
* Gj q's y “0D t 
XY vrs Ale re | 
aol ere 
y j I ? tly 
(arlene Pere 
Vy 
i 
iat rdt }- 
! be 2 F 


The exit sequence consists of checking X(58), (if set 
tie NONO table is printed), and KX(55), if set, (the exit 


clock time is printed). The NONO table codes for SUBE1 and 


SUBE2 are: 
O System can engage the target. 
1 The ASM failed at launch. 
2 fhe ship did not detect this target. 
3 The missile or raider was not identified as a 


target suitable for engagement. 
A the ship can not acquare the target. 
Reaction delay was too long. 


6 “The missi1lé will impact prior to-coming into 
max SAM range. 


Non missile ship. 
Detection will not occur until after end of game. 


Can not engage. 


10. SUBROUTINE SUBES 

This oF eee evaluates the engageability of the tar- 
Petean Light of Che current situation at the proposed launch 
fe,  oUBES is called by EFFECT, and it in turn calls 
meeon, STOREV, CPA, TIMMY, and, if X(55) is set, GETIMX on 
entry and exit. 

The inputs to SUBES3 are event time and a participant 
code. The code is an integer of the following type, IIIJJKLM, 
where III is the target number, JJ is the firing ship, K is 
the missile system, L the director, and M is the channel. 
number. The subroutine breaks down the code and then checks 
the appropriate portion of the SHIP table to determine 
whether the firing ship has been sunk, (if X(5) has been 


set, this check will indicate the ship has not been sunk 
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even if it has absorbed more hits than are considered nec- 
meaty totput it out oihvcaction)2 cIiiuthesship has “expended 
the magazine of this missile system, the subroutine exits. 
Having established the credibility of the firing ship, SUBE3 
determines whether the raider has been shot down or sunk, 
mada (SO, itcexits* 

SUBE3 now determines the maximum missile range for the 
system under consideration and compares this with the raider's 
CPA. If the target will not enter this maximum missile en- 
velope the subroutine exits. PRIOR is called to determine 
whether this ship/system/director combination has another 
shoot event at the present game time and, if so, if that tar- 
get is more threatening. PRIOR logic will consider the 
closest target as most threatening. If PRIOR returns the 
eeme target , the subroutine continues. * If not, it checks 
out the target as above and exits or continues as appro- 
priate. 

If the target has not been destroyed, SUBE3 checks the 
coordination set-up. If both the fleet and the firing ship 
are coordination capable, the subroutine checks the target 
availability. If the target is under engagement, the sub- 
routine calls STOREV to store a Type 3 event at the time the 
target will become available. If the target remains engage-— 
able, SUBE3 checks the NONO table. If the value O is found 
tie event continues. «If not; the subroutine exits. 

The subroutine now goes to the channel assignment phase. 


If the SAM or BLUE SSM requires no director, channel number 
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7 has already been assigned in SUBE1 or SUBE2, and the flag 
meOr is set«to 0... If the missile requires a director for 
panet or all of its flight, and’a channel has been assigned 
on a previous run through SUBE3, the ieroutide checks to 
see whether the assigned channel is busy. lied tays soLleol tires 
_ermco b> if not, ICOl is set to O. If the missile has not 
been assigned a channel previously, the availability of a 
channel is checked. A free channel, if found, is assigned 
and ICOT is set to 0; otherwise a search is made to deter- 
mine which channel will first become available and this one 
is assigned. If at this point all channels are busy and the 
director is required for the full SAM/SSM flight, STOREV is 
.called and a Type 3 event stored for the time the director 
will become available. 

rEtthe systems all indicate they are engaged in an 
intercept and upon searching the SAMASM .aand SAM tables no 
record of this can be found, the program terminates after 
printing an error message. 

SUBE3 now proceeds to check whether the firing ship 
itself is engaging the target. If it is the subroutine exits. 
If the target is an ASM or SSM, the subroutine. checks to see 
Eaatethe target will not impact prior to‘BLUE missile launch. 
me ec. will), the’ program, exits. 

At this point, if ICOT equals 1, the channels were all 
considered busy and the BLUE missile only requires the di- 
rector for terminal homing. SUBE3 checks the SAMASM listing 
to determine when the channel will be available and stores a 


Type 3 event. 
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Eeyicor rene 0, SUBE3 calls TIMMY to determine the 
intercept position, single shot. kill ee SAM/SSM 
fereht time, ‘and the intercept possibility given the BLUE 
missile flight and warhead characteristics. If TIMMY re- 
tuans the possibility flag -1, the intercept is not possible, 
but may be in the future. SUBE3 stores a Type 3 event for 
present time plus six seconds. If the indicator is zero, 
the intercept is not and will not be possible, and the sub- 
moubaine exits. If the indicator is one, the intercept is 
now possible and the program continues. If the intercept 
was possible, but the intercept point was in the hazard zone 
of another BLUE unit, TIMMY returns the indicator two. In 
this case, SUBE3 again stores a Type 3’ event for present 
time plus six seconds. 

If the SAM/SSM requires the director for only the ter- 
minal phase of flight, SUBE3 determines whether the director 
can be scheduled for the intercept. If not, the subroutine 
stores a Type 3 event for the time the director will be 
free, minus flight time, plus terminal homing mee plus 
0.06 seconds. 

If scheduling is not needed, or the director will be 
Milable when needed the missile is considered launched. 
The SAM-ASM-RAIDER and SAM tables are filled in and maissiz les 
on board and in the fleet counters are decremented. If 
eaere iS coordination, and the firing ship is coordination 
eapapblie, the target availability is up-dated. The launcher 


reload time is determined and stored. The channel 
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Bvatlability indicator is set to 1 and a Type 5 event is 
stored for intercept time. | 

Praor tor exact ,) Gifex(9) 1s set ,/ SUBES! prints out the 
firing ship number, the target number, and the reason for 
exiting the subroutine are printed. The breakdown for the 
REASON code for SUBE3 follows: 


OB bueSailvioi £ iired:. 

tf firing Ship sunk or Out of aciLion. 
2 Target being engaged. 

Target already destroyed. 


Firing ship out of missiles. 


oa ® Ww 


Director required for entire missile flight 
and is not available. 


6 Firing ship already engaging the target. 


7 Missile requires director only part time but 
all channels are busy. 


8 NONO table prohibits engagement. 

9 Launcher not available. 

10 ASM will impact prior to SAM launch. 
i1..,.Director cannot be scheduled for the intercept. 
12 Darector operation in horizon scan mode, busy. 


13 Target will not close BLUE missile envelope or 
TIMMY indicates no intercept. 


14 TIMMY indicates target may enter envelope later. 


ts LEVMY andicates tarset intercept point is too 
close to another BLUE unit. 


11. SUBROUTINE SUBE4 

SUBE4 determines whether the RED missile failed in 
meaeht, and if not, increments certain counters indicating 
Ship hits and possible sinking. SUBE4 is called by EFFECT 
aoa at in turn calls RANDOM, and, if X(55) is set, GETIMX 


en entry and exit. 
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SUBE4 first breaks down the input code into the RED 
ASM/SSM number, and the BLUE target ship number. The code 
enters SUBE4 as the integer RRRBB. ‘It then checks to see 
whether the ASM/SSM has been shot down, or had its launch 
negated due to the loss of the launching raider prior to the 
metual missile launch time. If neither is the case, the 
probability of successful flight is drawn from the X array 
and RANDOM is called to produce a random number. If it is 
peeeininied that a failure occurred, the ASM table is set to 
3 and the subroutine exits. Otherwise the target ship is 
assumed to have been hit and the appropriate counters are 
| incremented. If X(5) is not set, the subroutine compares 
the number of hits taken to the amount of damage the ship 
can withstand, and if too many hits have been taken, de- 


clares the ship sunk. 


12. SUBROUTINE SUBES5 

SUBES evaluates the BLUE missile intercept and performs 
checks on other ASM/SSM if the target shot down or sunk was 
are raider. ‘SUBES is’ called ‘by EFFECT and it im turn 
eabis STOREV, REACT, RANDOM, and if X(55) is set, GETIMX on 
entry and exit. 

SUBES is called with an input variable of salvo number. 
Using this, SUBES5 gets the target number, firing ship number, 
missile system, dit Soo RAIAnd channel number from the SAM- 
ASM-RAID table. The director reassignment and reacquisition 


delays are determined and stored, fleet and individual ship 
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statistical datatare taken, and a check is made to insure 
that he target is still engageable. If the target has al- 
ready been sunk or shot down, the “SAM-ASM-RAID table is set 
foe—-1, indicating an overkill. If the target is an ASM/SSM 
phat failed in flight or has already impacted, the salvo is 
determined to have been a failure and the SAM-ASM-RAID table 
He Set to zero. If the target is engaged, the probability 
Sf£eBLUE missile successful flight and support is pulled 

from the X array, and the single shot kill probability is 
removed from the SAM-ASM-RAID table. RANDOM is called to 
provide a random number. Success is determined by comparing 
PK to the random number drawn. This is done by determining 


Pre as follows: 


PE =_P 
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Where PS is the probability of successful flight, PKSS 
is the single shot kill probability, and N is the number of 
BLUE missiles in the salvo. If the random number is less 
than PK, the shot was:‘successful. If the shot was unsuccess- 
ful and the missile requires a director, the channel is re- 
leased and Type 3, shoot, event is stored for present game 
bame. If the missile doesn't require a director, just the 
Shoot event is stored, and the subroutine exits. 

If the shot was successful and the target was an. ASM/ 
SSM, the SAM, SAM-ASM-RAID, ASM, and kill counters are up- 
dated, a Type 6, channel free, event is stored, and the 
Eubroutine exits. If the target destroyed was an aircraft, 


ae at the number of hits has been sufficient to disable a 


116 


; 7 
t 


Shunt! RAND 


£0 Seems sk 1 ~~ 5 foe ee 


ee : f it web soe vs 


vo waea- 





surface raider, a following stevs are taken prior to the 
dar table update and exit. A check is made to determine 
maether the raider has any ASM's or SSM’s aboard. i1f so, 
the ASM counter is incremented and the appropriate numbers 
are entered in the ASM table indicating that these missiles 
were lost with the launching platform. Next the ASM table 
is checked to locate any missiles launched by the raider 
just destroyed. When one is found, its launch time is com- 
bared, tO raider kill time, and if the launch time is greater, 
entries are changed in the ASM table to correspond. Upon 
completing the search, the program continues to exit where 


mamecnecks X(55),,X(56).,. and. X(58),.for debug. output. 


13. SUBROUTINE SUBE6 

SUBEG releases the fire control radar channel at the 
time of the event. SUBE6 calls GETIMX upon entry and exit 
He m(oo) is set. 

The input value to SUBE6 is an integer code AABCD, 
where AA is the firing ship number, B is the missile system, 
C is the director number, and D the channel number. After 
breaking down the code, SUBE6 sets the appropriate portion 


of the CHANEL array to zero, and exits. 


14. SUBROUTINE TIMMY 

TIMMY computes the predicted intercept point and flight 
time for the BLUE SAM/SSM and determines whether the inter- 
cept lies within the intercept contour. TIMMY is called by 
SUBE3 and TIMMY in turn calls TMDASM, ASMTIM, RADIAN, ZLNKA, 


mOCTIM, and if X(55) is set, GETIMX on entry and exit. 
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Inputs to TIMMY are target number, firing ship number, 
missile System number, present target position, and present 
game time. Outputs are single shot kill probability, flight 
time, and an indicator flag which says the “nl eroent is pos- 
Sapte, is not and will not be possible, or-is not but may be 

spoessible in the future. 

Using the target number, TIMMY decides whether the tar- 
get is a raider or an ASM/SSM and determines the end of game 
position for raiders and target ship position for an ASM/SSM. 
If X(8) is set, firing sectors have been assigned and TIMMY 
checks to ensure that the assigned firing sector and the 
faunehner train limits permit.a launch: If launch is pro- 

eerprted, REASON is set to 2, flag IND is set:to zero, and 


the subroutine exits. li Ghe® dauneh 


there ma 


is possible ay 


? 


be aS many as two firing zones. These two occurrences are 


Shown in Figure Bl. 
SECTOR 
LIMITS 















LAUNCHER TRAIN LIMITS LIMITS 


LIMITS 
LAUNCHER 


Figure Bl 


One or Two Firing Zones 
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Using the negtlEs of the ZLNKA reordering of the sector 
mmrcs and launcher train limits, TIMMY determines whether 
there are two zones or one, and which of the bearings con- 
Pmurute the actual firing zone, or zones, Jamits: If no 
Sectors have been assigned to the fleet, or this unit, the 
firing zone limits coincide with the launcher train limita-— 
tions. 

TIMMY next checks target :-altitude against the maximum 
SAM altitude. lt the target isitoo high, REASON’ is set to 
eleven and the subroutine exits. Otherwise, the intercept 
parameters are determined. These intercept parameters are 
the values entered into the X-array that approximate the 
actual intercept contour. The intercept parameters are A, 
meee SLOPE, and SEEK. The actual intercept contour depends 
on target size, speed, and altitude as well as missile pa- 
rameters. There are usually more than one contour for any 
given missile system. The actual contours are found in the 
various technical publications dealing with the specific 
missile. OP-3594, Performance and Installation Characters 
istics Surface Missile Weapon Systems, also contains general 
contour nomograms, and the Applied Physics Laboratory of 
Johns Hopkins University has developed computer simulations 
which generate contours. A typical contour, along with the 
TIMMY approximation, is shown in Ficunme BBR A is the off= 
range capability, B is the down-range capability, and N is a 


degree of curvature. 
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Figure B2a 
The Actual Intercept Contour 
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Figure B2b 


The Approximation of the Contour 
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The maximum missile range for the given target altitude 
is determined, and if the target will not enter this range, 
as determined by TMDASM, TIMMY exits after setting REASON to 
seven. If the target is an ASM/SSM and the target ship is 
ime firing ship, target course is found and SAM flight time 
is determined by: 


FT = DLT/(SSPD+VEL) 


where DLT is the distance between the firing ship and the 
target, SSPD is the speed of the SAM, and VEL is the speed 
ef the ay ae nen In this case the PK that is returned to 
SUBE3 is that for a non-crossing target. 

i the tarset ship asinot) the: firing ships), or the) tar- 
get is a raider, the following angles are determined: ALPHA, 
the angle between the target flight path and the line o£ 
Sight at launch, and THETA, the angle in Cartesian coordi- 
hates of the present line of sight. The quantity GSINE is 
determined by: 


GSINE = SIN(ALPHA)*(VEL/SSPD) 


If the absolute value of GSINE is greater than one, the 
intercept is impossible and TIMMY exits after setting REASON 
equal to ten. Otherwise, the angles GAMMA, BETA, and DELTA 
are computed. 


GAMMA 


ARCSIN(GSINE) 


BETA = 180° - GAMMA - ALPHA 


DELTA 180° - BETA 
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GAMMA is the angle between the SAM/SSM flight path and the 
line of sight, BETA is the angle between the SAM/SSM flight 
paeh and the ‘target flight path, and DELTA is the compliment 
Sipe rA., if XC59) is) set, “the values ‘for ‘target: number ; 
firing ship number, ALPHA, BETA, GAMMA, present target posi- 
tion, GSINE, and present range to the target are printed. 

PK for this type of target is determined by DELTA. If DELTA 
moepreater than 45°, the target is declared a crossing tar- 
get and the appropriate PK is retrieved from the X-array. 
Otherwise, the probability of kill for a non-crossing target 


is retrieved. Flight time is determined by: 
FT = DET * SIN(ALPHA)/(SSPD * SIN(BETA) ) 


Using present time plus flight time for the input, TIMMY now 
calls LOCTIM or ASMTIM, as appropriate, to determine the 
intercept coordinates. The range to the intercept is com- 
puted and compared to the minimum L MeeiTS range. If the 
range is too short, IND is set to minus one, REASON to four, 
and the subroutine exits. The launch bearing is calculated 
and the first of two possible sets of firing zone bearings 
are compared to determine whether the launch bearing is with- 
mwaetune firing bearings. If it is not, TIMMY checks to see 
if there was a second firing zone, if so it checks this one, 
otherwise, IND is set to minus one, REASON to three, and 
TIMMY exits. If the intercept angle is good, the intercept 
parameter is checked. A, B, and N are combined to form the 
equation: 


oi GE) = 
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Pevalue of twontor,N results: inea icircleor an ellipse. A 
walue of six resulisiin a near rectangle. Fractional values 
for N are permitted. 

The "ellipse" is translated so that fee center is on 
moe firing ship, and rotated so that the positive Y axis is 
apposite to the line of flight of the target. The intercept 
coordinates are substituted into the equation, in absolute 
value, and if the resulting number is less than or equal to 
one, the intercept will be within the ellipse. The next 
check is the seeker limit test. The y intercept value and 
the absolute value of the x coordinate are entered into the 
equation: 


Z = YI — SLOPE * XI 


If Z is greater than or equal to SEEK, the problem continues. 
If not, TIMMY exits after setting REASON to six. The sub- 
routine now checks the X-vector to determine whether there 
are any restrictions on how close an intercept can be to 
another friendly unit. If there are, the hazard zone is 
eonstructed. If the intercept falls within the hazard zone, 
REASON is set to nine, IND is set to two, and TIMMY exits. 

If the problem has progressed this far, the intercept 
is considered possible and the BLUE missile salvo is launched 
while IND is set to one. If the target is an ASM/SSM and 
mae firing ship is the target, the program exits. If not, 
the target altitude is figured into the flight time prior to 
feucurn to SUBES. If X(9) is set, TIMMY prints out the ship 


number, target number, IND value, and REASON prior to exit. 
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Reason codes as used in TIMMY 


ae Not used. 

2 No firing zones for this system. 

3 Firing bearing outside sector. 

4 Lareet too close. 

5 Target intercept point outside engagement contour. 
6 Target beyond seeker limits. 

i Target beyond range and outbound. 

8 Intercept after ASM/SSM impact. 

9 Target too close: to BLUE unit at, antercept. 
10 GSINE greater than one. 

if Target too high. 


15. SUBROUTINE JAMM 

The subroutine JAMM determines the range at which a 
monostatic radar can detect a given target in a multiple 
standoff and self-screening jamming environment. 

The calling routine, either SUBE1 in the case of raiders, 
or SUBEH2 in the case of missiles, first. determines that 
there are jammers in the problem and that at least one jam- 
mer is on the same frequency band as the detecting radar. 

If both situations obtain, subroutine JAMM is called. 

JAMM initially checks to see if the target will close 
the detecting unit within the minimum of the clear environ- 
Ment, radar horizon, and theoretical range of the radar. If 
the target will not close within this minimum range, the 
Subroutine is exited with a flag set to show that detection 
would not occur. In the case where the target will enter 
this minimum range, JAMM computes the time at which the tar- 


get will arrive at the closest point of approach (CPA) to 
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the detecting radar. Appropriate checks are made for sta- 
LPonary targets and targets with only minor changes in their 
mor Y position components. JAMM next obtains the bearing 
from the radar to the target at the CPA BsoaP ibs for use in 
selecting the appropriate side-lobe coefficient. Individual 
radars are then checked to see if they have been shot down 
perer to the target arriving at CPA (applicable only for 
missile detections), if they are jamming on the frequency 
band of the detecting radar, and if they are above the de- 
meceing radar's horizon. if -a2illorrihe conditions are met? 
the bearing from the jammer to the jammer's jamming target 
is found in order to determine whether the detecting radar 


is within the main beam or the side-lobe of the jammer. 


jammer is then computed, taking into consideration the power 
of the jammer, the range from the jammer to the detecting 
radar, and the side-lobe/main beam relationship between the 
detecting radar and the individual jammer. The jamming power 
arriving at the detecting radar is then summed over all of 
the jammers and used to compute the burn-through range for 
ome radar against the specified target. Prior to exiting 
Ehbe Subroutine, the range to the target's CPA to the detec- 
tines radar is checked to ensure that it is less than the. 
burn-through range. If the distance to CPA is greater than 
the burn-through range a flag is set indicating that detec- 
tion did not occur due to CPA considerations. Otherwise the 
Subroutine returns as the detection range the minimum of the 


elear, jamming, radar horizon, and theoretical ranges. 
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16. SUBROUTINE BIDET 

This subroutine computes the number of detection zones, 
ranees, and times of detection for a bistatic radar in 
either a clear or jamming environment. BIDET is called by 
SUBE1 and SUBE2, and it in turn calls GETIMX on entry and 
exer iid) X(55) is set. 

BIDET begins by setting the initial values of aircraft/ 
pow start position, final position, altitude, and velocity. 
The transmitter and receiver positions and. radar antenna 
height are established for the bistatic system under consid- 
eration. A check for non-moving targets is conducted prior 
to entering the closest point of approach (CPA) routine. In 
. the case where the target has moved less than 0.1 nautical 


mae an both the X. and Y direction 


, the target's speed is 
considered to be Zero and the target enters a special loop 
to determine whether it is detectable at its initial posi- 
tion within the constraints of both radar horizon and radar 
Signal strength. 

For the moving target, the CPA of the target's projected 
track to the transmitter is determined by solving simultane- 
ously the equations of the target's track and the line per- 
pendicular to the target's track and passing through the 
transmitter location. The CPA to the receiver is solved in 
a similar manner. If movement in either the X or Y direc- 
tions is less than 0.1 nautical mile, the CPA is established 
py setting’ the CPA coordinate for the small changing axis to 


the target's original position corresponding coordinate and 
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the alternate coordinate to the corresponding alternate co- 
ordinate of the transmitter or receiver. _ For example, if 
the X coordinate of the target's positions changes less than 
0.1 nautical mile, the CPA to the receiver is set to the X 
coordinate of the target's initial position and the Y co- 
ordinate of the receiver's position. 

hhe Gadar horizon for the, transmitter and the receiver 
are then compared to the distance to their respective target 
track CPAs to determine whether the target enters the radar 
horizon of both the transmitter and the receiver. If it 
fails to enter both, the routine is exited with a no detec- 
tion. Otherwise, the times that the target first appears 
above and disappears below the transmitter and receiver ra- 


dar h 


orizons are computed. The time interval that the tar- 
get is simultaneously above both horizons is determined by 
taking the maximum of the appearance times and the minimum 
of the disappearance times. If the target is never simul- 
taneously above both horizons; or the appearance time is 
greater than maximum game time; or the disappearance time is 
prior to start of game time, the subroutine is exited with a 
no detection. 

Next, detection times for the target based on radar 
Signal strength are computed, and the intersection of these 
time intervals with the above joint radar horizon time in- 
terval determines the time during which the target is de- 
tected. No detection occurs if the intersection of the two 


time sets (i.e. radar strength detection and radar horizon 


e@ervection) is the null set. 
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If a detection has occurred in the clear environment, 
the subroutine checks to see whether any jammers are in use 
and if they are on the same frequency band as the bistatic 
radar. if jammers are in use, the time that the target 
errives at its midpoint in the bistatic detection envelope 
is determined. This time is used in arriving at an average 
jammer position to compute the jamming power at the radar 
during the detection phase. This will tend to provide a 
somewhat conservative estimate of detection ranges in the 
face of jamming because the jammers are in general going to 
be nearer to the radar and the range square losses of the 
jamming power will be reduced. In the case of the station- 
ary target, midpoint time is set to an arbitrary three min- 
utes after start of game time to allow for some movement cf 
the jammers tomar the force. The subroutine then adds the 
jamming power from each jammer... It does this by checking 
all of the raiders to determine if the individual raider has 
been shot down prior to the midpoint time, and if the jammer 
is above the radar horizon and is jamming on the frequency 
band of the radar. If these conditions are met, the sub- 
routine examines the geometry of the jammer's position and 
the radar to determine the relationship between the jamming 
and radar antennas. With this relationship, the appropriate 
Side-lobe/main beam antenna coefficients are used to adjust 
the jamming level. 

The constant in the bistatic range equation is then ad- 


justed to reflect the jamming environment and the computation 
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previously conducted in the clear environment are repeated 
to determine whether detection will occur in the jamming 
situation. 

The subroutine exits with the number of detection zones 
woich occurred; the time and range of initial detection in’ 
wae first zone, the time of last detection in the first 
zone; and, if applicable, the time and range of first detec- 
tion in the second zone along with the time contact is lost 


in the second zone. 


17. SUBROUTINE DTIMD 
Subroutine DTIMD is a double precision subroutine that 

determines the real and imaginary roots of the fourth order 
polynomial which is used to solve for the bistatic radar de- 
tection times. The subroutine employs the Newton-Raphson 
iterative technique with the final iteration on each root 
performed using the original polenenia’l rather than the re- 
duced polynomial to avoid accumulated errors a the reduced 


polynomial. The subroutine was obtained from the Naval 


Postgraduate School Computer Center library. 


18. SUBROUTINE CPA 

Subroutine CPA obtains the closest point of approach of 
a target's track to a reference point by solving the simul- 
taneous equations of the line representing the target's 
track and the line through the reference point and perpen- 
Gieular to ther target's track); Checks are conducted for 


target tracks which pass through the reference point, 
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represent stationary targets, or targets which have only 


minor changes in either the X or the Y direction. 


19. SUBROUTINE REACT 

Subroutine REACT allows the use of reaction time delays 
marouch the use of probability distribution functions. This 
subroutine iS Similar to, if not identical to, one developed 
for use in FLOATS by APL at the Johns-Hopkins University. 
The distribution functions may be entered for seven classes 
et edelay. They are: 


1 Target detectability to actual detection. 


2 Detection to ID and designation to the weapons 
direction equipment. 

3 Director assignment delay. 

4 Director acquisition delay. 


THs a+ 
4 


ijrane delay 
di Lis iwiay. 


6 Target detectability to detection for horizon 
Scan mode. 


i Decision delay in horizon’scan mode. 

Within these class areas, up to five separate types may 
be designated. REACT uses the type and class numbers to de- 
termine a delay time which it returns to the calling sub- 
routine. REACT is called by SUBE1, SUBE2, and SUBES. 

REACT in turn calls only RANDOM. 

The inputs to REACT are the class number and type, and 
the output is the delay in minutes. The class number is. 
used to locate that portion of the X-array where the delay 
is stored. The type number is used to localize the address 
within the area. In employing its Monte Carlo technique, 


REACT calls RANDOM and converts the returned uniform (0,1) 
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number to an aan mi Lf this sntecer as 10, the subroutine 
Secs it equal to 9. If the integer is zero, the RANGE, de- 
fined as the difference in time between two adjacent blocks, 
ms set equal to the value in the first block. Otherwise, 
the address is determined using the class and type numbers 
and then adding on the integer. The delay consists of two 
wees. Lhe first part-is the time in the address biock. 

The second part is determined by drawing another random num- 
ber and multiplying it by RANGE. Except in the case where 
the integer is zero, RANGE is determined by subtracting the 
time in the address block from that in the next highest 
Brock. in the “zero case, the total delay consists of this 
eeaection Of RANGE. In all other cases,.the fraction of 
RANGE is added to the time stored in the address block to 
Obtain the total delay. At this point REACT returns to the 
Salting routine. 

As an example of a distribution function, consider the 
meicine time for an aircrait for a catapult launch. Data is 
available for one hundred aircraft. The maximum Gahcine 
time was seven minutes. The distribution of waiting times 


was as follows: 





jength of wait Ql) ea a8) ae ee eo 
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Figure B3 


An Example of Waiting Time 
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The plot of length of time waiting against number of occur- 


rences would look like this: 
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Figure B4 


Histogram of Waiting Times 


The probability of waiting between 0.0 and 1.0 minutes is 
0.1, the probability of waiting between 1.0 and 2.0 minutes 
o0.4. When the histogram is transformed into a cumulative 


probability distribution function, the following results. 
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Cumulative Distribution Function of the Waiting Times 
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This distribution is what is entered in the X-array. 
The first block would contain the time for .1 cumulative 
probability. This would be one minute. Block two that for 
femeumulative probability, or 1.25 minutes, and so on up to 


block ten which would have entered a seven for seven minutes. 


20. SUBROUTINE PRIOR 

Subroutine PRIOR is called by SUBE3. PRIOR calls REMOVE, 
SIOREVY, LOCTIM, ASMTIM, TARGET, and if K(55) is set, GETIMX 
on entry and exit. The purpose of PRIOR is to determine the 
most threatening target among those that a system may choose 
from at present game time. 

PRIOR first computes the value CHECK for the target 
being compared. CHECK is defined as maximum missile system 
range divided by present target range. ‘TARGET is called to 
provide PRIOR with a candidate for comparison. The same 
measure of threat is obtained for the candidate target and 
Stored in CHERRY. If the present value in CHECK is greater 
than or equal to CHERRY, PRIOR goes to the next candidate, 
mebnwere is one. If CHERRY is greater than CHECK, the val- 
ues of target number, channel number, and target position 
are saved. CHECK is made equal to CHERRY and PRIOR seeks a 
new comparison. 

If the value returned by TARGET in MORE is not zero, 
there ie another candidate to be evaluated. . If MORE is zero, 
there are no more comparisons to make and the final section 
of PRIOR is entered. The stored value of most threatening 


target number is compared with the target sent to PRIOR 
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Srafinally from SUBE3. If they are the same, the subroutine 
exits. If PRIOR has found a more threatening target, REMOVE 
is called to take the type 3 event associated with the new 
harget out of the events list.  STOREV is then called to 
enter a type 3 event for the target being intercepted. 

The debug switch X(9) is checked at exit and if it is 
set, old target, new target, missile system (in the code 
AABC were AA is ship number, B is system number, and C is 
director number), old channel, and new channel are printed. 


X(55) is checked and PRIOR returns to SUBES. 


21. SUBROUTINE RCHECK 

Subroutine RCHECK evaluates the prospective target ship 
to determine if the raider can in fact "see" the target. 
REHECK is called by SUBE2 and RCHECK calls GETIMX on entry 
pme exit if X(55) is set. The inputs to RCHECK are the 
raider number and proposed target ship number. The two di- 
mensional distance between the two units is computed. The 
antenna heights on the target ship are compared and the 
tallest one is used to determine the ship's radar horizon. 
If the ship has no radar, a height of 15 feet is entered. 
The combined radar horizons, maximum theoretical radar range 
of the raider's radar, and the actual range are compared. 
it the target is within both the radar horizon and the maxi- 
mum one range, the go/nogo flag is set to zero. If either 


mest fails, the flag is set to one and the subroutine exits. 
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22. SUBROUTINE TMDASM 

This subroutine determines if the target, or other unit 
in question, will enter within a certain input range of a 
iweeunit;, and if so,.at what time will it enter and leave. 
TMDASM is called by SUBE1, SUBE2, and TIMMY. TMDASM in turn 
ealis GETIMX on entry and exit if X(55) is set. 

TMDASM begins by determining the range from the refer- 
Enee ship to the target unit. If this range is less than or 
equal to the input range, the time of entry is considered 
[aeperthe input time. If the range to the target unit is 
greater than the input range, TMDASM determines the slope of 
a line representing target unit movement. Using this and 
‘the present target unit position as a reference point, TMDASM 
forms tie iine equation. This equation and the equation for 
a cireie with ea ate equal to the input range are solved 
Simultaneously to determine the points of intersection. 
This is done by substituting the line equation into the 
circle equation and solving the quadratic equation. If the 
roots are imaginary or equal, the entry is either not made 
or the length of stay is infinitely small. These cases are 
hmegected as no solution. Upon finding. two real roots, ehe 
Subroutine determines which of the two points the target 
unit would reach first, and which one it would reach second. 
fmmene target) unit.is already further down its track than 
the exit point, no entry is indicated. Otherwise the time 
of leaving the circle is determined and, if the target unit 


was not already inside the circle, the time of entering. 
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23. SUBROUTINE TMDBOM 

This subroutine determines when the raider in question 
will reach a given range from the force center, position 
(0,0). TMDBOM is called by MAIN, and JULIE... TMDBOM in turn 
Gees) GETIMX on entry and exit if X(55) is set. 

TMDBOM proceeds in the same way as TMDASM except that 
the reference position, center of the range circle, is al- 
ways force center instead of a firing ship. TMDBOM only 
computes the time that the raider will cross a given range. 


It neglects the departure time. 


24. SUBROUTINE ZLNKA 

Subroutine ZLNKA is a copy of the subroutine SHSORT, 
mentioned below, as it existed on 1 February 1975. ZLNKA 
inputs a vector of-real numbers and an ordered vector of 
integers. ZLNKA reorders the real vector from smallest to 
largest and reorders the integer vector to the same order as 


the real vector. ZLNKA is called by TIMMY. 


29. SUBROUTINE TJSORT 

Subroutine TJSORT is a modification of the subroutine 
SHSORT which is available from the Naval Postgraduate School 
Computer library. .TJSORT is called by SUBE2, and BIDET. 
This subroutine inputs a vector of real numbers and rear-. 


ranges them from smallest to largest. 
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26. SUBROUTINE LOCTIM 

LOCTIM computes ‘the raider’ position: at time T. It is 
called by MAIN, SUBE1, SUBE2, SUBES, JAMM, BIDET, TIMMY, and 
PeeOok. LOCTIM in turn calls GETIMX' on entry and exit if 
hoop) is set. 

LOCTIM first determines raider velocity and the time of 
the raider position stored in the RAID table. The distance 
the raider travels between the present time and the stored 
time (DL) is computed. If the velocity is zero, present tar- 
get position is set to the initial: position and the program 
exits. The distance between the stored position and the end 
of game position (D) is determined and the ratio DL/D is com- 
puted. The X,Y position at the time desired is determined 
by: 


xPo= XZ = AX +t DI 


XP is the X position at the time of interest 
YP as the Y position at the time of interest 
XZ is the stored X position 
YZ is the stored Y position 


AX is the distance between XZ and the end of game 
position 


AY is the distance between YZ and the end of game 
position 


DE as the vatio Diy D 
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If the input flag IND has been set.to 5, the new time 
and position are stored in the RAID table, otherwise the 


subroutine exits. 


27. SUBROUTINE ASMTIM 
ASMTIM computes the ASM/SSM position at time T. It is 
called by SUBE3, TIMMY, and PRIOR. ASMTIM calls GETIMX on 
eaceyiand) exit) 1f)X(55) asset. 
The toga ceof, ASMTIM is adentical to that. of LOCTIM but 
the definition of variables has been changed. 
XZ becomes the ASM/SSM X launch coordinate 
YZ becomes the ASM/SSM Y launch coordinate 


AX becomes the distance between the launch position 
and the X coordinate of the target ship 


AY becomes the distance between the launch nosition 


and the Y coordinate of the target ship 
The ASM/SSM position is now computed in the same manner 
as in LOCTIM and the subroutine exits. The positions are 


not saved by ASMTIM. 


28. SUBROUTINE RANDOM 

Subroutine RANDOM is a generator of Uniform (0.0, IQ) 
random numbers. It was programmed by Gerald P. Learmouth, 
of the W.R. Church Computer Center, Naval Postgraduate 
penool;, The following is taken from the library description 


of the subroutine: 
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"The basic generator is a Lehmer congruential 
eof the form: ; 
Che) = AACN ye ieDr CP") 

where A = 7**5 = 16807 and P = (2**31)-1. The 

resulting X(N)'s are uniformly distributed random 

integers between 1 and P. The uniform REAL*4 num- 

bers are formed by right-shifting the integer X(N) 

and appending a proper floating point exponent to 

form a number between 0O.O and 1.0." 
29. SUBROUTINE HISTG 

Subroutine HISTG is found in the subroutine library of 
the Naval Postgraduate School subroutine library. A copy of 
HISTG has been included at the end of the program. The sub- 
routine is used if the user wants a statistical display of 
the number of hits a certain specified BLUE unit will take. 
Also provided is a statistical display of the time interval 
over which the RED missiles arrive, and one of the hits per 
Minute. A histogram of the selected parameter is provided in 
all cases. HISTG is called by MAIN just prior to run com- 


pletion if X(85) has been set and the conditions detailed in 


section V.B have been met. 


30. SUBROUTINE RAIDIS 

This subroutine provides a plot of the raider tracks 
mom their initial position to launch position. The maximum 
Saze of the plot is 800 nm from force center. In the case 
where a raider does not launch weapons, the 30 minute game 
time position is used ag the terminal point of the track. 
woethe launch or 30 minute position is more than 800 nm from 
moree center, a label of “no launch" is displayed at the 


Baraer's initial position. 
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oa. SUBROUTINE DISPL 

DISPL is similar to RAIDIS except that the BLUE unit 
positions are plotted and they must’ be within 200 nm of 
force center. A 100 nm range circle is included in the dis- 


play. 


32. SUBROUTINE SETIMX 
inas Subroutine sets the internal clock to zero. The 


entry GETIMX gives the internal clock time when called. 


33. FUNCTION RADIAN 
The function RADIAN converts degrees to radians. It is 
called by MAIN, JULIE, and TIMMY. The conversion is done 
with: 
RADIAN = RD * 0.017453293 
where RADIAN is the output angle in radians and RD is the 


input angle in degrees. 


140 


Pe rar Ol we rl itor a ' 





APPENDIX C 


RADAR DETECTION RANGE INPUTS 


This appendix provides the equations used to determine 


the 8 and a inputs for the individual radar classes. 


1. CLEAR ENVIRONMENT 
B is the range at which a one Square meter target can be 

detected in a non-jamming environment. Kets 2 constan’ Ot 
the radar system and in the monostatic case provides the 
range of detection for the radar by the relation: 
Range = B(o)? where o is the target's radar cross section. 
“The general equation for 8 is: 

Pee 4 

Be ib Yea, 58! 
m 





where Pi is the power of the eeaMenaa ber a watts; G is the 
antenna gain coefficient; »X is the radar wave length in 
meters; and a is the minimum detectable signal of the radar 
with appropriate consideration for probability of detection 
and propagation losses. If the above units are used in 
calculating 8, the result will be a detection range in meters 
which must be converted to nautical miles for use in the 
Simulation. Values of 8, in nautical miles for a one square 


meter target, are reported in Ref. 9 for standard U.S. Navy 


radar systems. 
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in the case of the bistatic system, .the equation for 6 


becomes 


where the possibility of different gains of the transmitter 
and receiver must be considered and the propagation losses 
in aan must take into consideration the differences in the 
transmitter to target path and the target to receiver path. 
Badationally, in the bistatic situation, the detection 

range is non only a function of the range from the receiver 
(R,.) To the? targetebue at is also a factor of the range from 
the transmitter (Ri) to the target. Thus the bistatic range 


equation becomes 


The monostatic radar cross section provides a suitable esti- 


maue for the bistatic radar cross section [Ref. 10, p. 592]. 


2. JAMMING ENVIRONMENT 

a is the range at which a one square meter target can be 
detected in a jamming level of one watt/mhz. The burn through 
range for a monostatic radar in a multiple-standoff jamming 


environment is: 


1 1 
a bee a 
Rds to>(oyiSUP) 


where o is the target's radar cross section and SJP is the 
sum of the jamming power arriving at the radar. The jamming 


power from an individual jammer on the radar's frequency is 
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a function of the power of the jammer; the side-lobe/main 
beam relationships between the jammer and the radar; and the 
square of the range from the jammer to the radar. The equa- 


tron is: 


Sue ie ID ee 
i qe at a 


where C; is the side-lobe coefficient for jammer i and the 
radar ; JP, is the jamming power of jammer i on the radar's 
frequency; and R; is the range from the jammer to the radar. 


a is found by the relation: 


where P, is the power of the transmitter; G is the gain of 
the radar (In bistatic systems G is the square root of the 
product of the Panini: and receiver gains. ); By LS: the 
meeeiver bandwidth of the radar; and S/J is the signal to 
jamming ratio for the desired probability of detection and 
appropriate propagation losses. In the bistatic range equa- 
tion in a multiple standoff jamming environment, the burn 


through range is found in a manner similar to the clear en- 


vironment to be 


au 
Pen = nace SIP). 


where Rit is the burn through range, Ri is the range from 
the bistatic transmitter to the target and the remaining 


mectOors are: as’ defined above. 
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APPENDIX D 


PROGRAM ARRAY LISTING 


This appendix provides an alphabetical listing of the 


principal arrays used in the computer program. The listing 


contains the values of the arrays for the indicated array 


arguments. The arguments for the arrays are specified in 


the heading of their respective listing. The dimensions of 


the individual arrays are stated in the COMMON statement at 


the beginning of each subroutine. 


ASMT (I,J) 


The ASMT array contains data concerning the RED missiles. 


I AMS/SSM launch seauvence number 


J=1 ASM/SSM status code as follows: 


0 


abkwnre 


eu Th 


Launch failure 

Successful launch 

Killed by SAM 

Inflight failure 

Hit target ship 

Launch platform shot down or sunk 
prvor’ to Launch 


J=2 Time of impact on target ship 


J=3 Raid number of launch platform 


J=4 Target ship number 


J=5 Time of launch 


J=6 Weapon class 


J=7 X coordinate of launch position 


J=8 Y coordinate of launch position 
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CHANNING gd (Ke Ls) 


An array used to keep track of which channels are busy. 


i Ship number 
J System number 
K Director number 
L Channel number 
O = Free 
1 = Busy 
JACK (CIES) 


This array stores the engagement and detection statistics 
for individual ships. 
I Ship number 
J=1 Chances to engage. The number of targets that 
will enter within the maximum missile range of 
the ship 
=2 Number of engagements 


J 
J=3 Number of targets destroyed 
J 


—4 ‘The number of raiders’that will enter the mis- 
Sile envelope. 


J=5 The number of raiders that will enter the en- 
velope and are detected 


J=6 The number of ASM/SSM that will enter the 
envelope 


J=7 The number of ASM/SSM that will enter the 
envelope and are detected 
NADAR (I,J) 
This array indicates the operational status of the on- 
board BLUE equipment. 


I Ship number 
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Search radar number 1 


O = down 
1 = up 


search radar number 2 


Search radar number 3 


Battery one, fire control radar number 1 


down 

operational 
no director 
operational 


normal mode 
needed for missile 
in horizon scan mode 


WNHO 


Battery 
Battery 
Battery 
Battery 
Battery 
Battery 
Battery 
Battery 
Battery 
Battery 


Battery 


one, fire 


one, fire 
twos’ fire 
two.” Lire 
two? fre 
three, 
three, 
three, 
OUI, 


POUR, 


TOU, tre Control radar 3 


control 


Control 


control 


Control 


contro! 


Launcher for battery one 


0 
a 


oll 


non operational 
operational 


Launcher for battery two 


Launcher for battery three 


Launcher for battery four 
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RAID "C1 JS) 
This array stores information on the raider's location 
and weapons. 
I Raider number 
J=1 me viector starting index of “the raid -ellass 


J=2 Time that raid was at position indicated in 
J=3,4 


J=3 KX ‘coordinate at tame T 
J=4 Yecoordinate “at time T 
X coordinate at end of game 
J=6 Y coordinate at end of game 
J=7 X coordinate at launch if the raider has weapons 
and can reach the launch position, otherwise 
the X coordinate at end of game. 
J=8 Y coordinate at launch, same as J=7 
J=9 Raider altitude in feet 
J=10 Raider's radar horizon 
J=11 Weapon class of raider's longest range weapon 
J=12 Total number of weapons assigned 
J=13 Status 
O = destroyed or sunk 


number of hits the raider may take until 
destroyed 


Z 
| 


J=14 X coordinate of raider at end of game or X 
coordinate of target ship for ASM/SSM launched 
by taider I 

J=15 Y coordinate comparable to J=14 


J=16 Time raider will reach second launch point if 
it is carrying two types of missiles 


J=17 Time the raider was destroyed 
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SAMASM (I,J) 
This array contains information covering the BLUE missile 
launch, target, engagement outcome, and provabidaty of kati. 
I BLUE salvo number 


J=1 Target number, 1-90 are raiders, 91-191 are 
ASM/SSM 


J=2 Ship number followed by system number 

J=3 Director number followed by channel number 

J=4 Time of launch 

J=5 Engagement outcome. Negative number indicates 
an overkill, zero a failure, and a positive 
number a successful shot. 


J=6 Time the channel will be free for reassignment 


d-7 YSangle shot probabihity .o£ kill 


SHIP (1i,J) 
The SHIP array contains information concerning the BLUE 
units. 
I Ship number 
J=1 <Ship,.cilass 
d=2 sxX,coordinate,of+ship, position 
dJ—3 o¥.coeordinate of» Ship.position 


J=4 Number of hits ship can receive until declared 
OUbsOnweriton. li Oo. Shap as out tol action. 


J=5 Ship class starting index 
J=6 Missile status code 


0 no missile capability 
if missile capability 


tl 
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SMAT (I,J) 
The SMAT array lists the results of the SAM firings. 
I SAM number 
J=1 SAM class index 
J=2 Firing ship number 
J 


=3 Number of missiles remaining in the battery 
firing the SAM 


J=4 lime Lire control radar channel 1 free 
J=5 Time fire control radar channel 2 free 
J=6 Time fire control radar channel 3 free 
J=7 Time fire control radar channel 4 free 
J=8 Time fire control radar channel 5 free 
J=9 Time fire control radar channel 6 free 


J=10 Number of SAMs in the salvo 


SATS Gh as) 
The STATS array is used to compile statistics concerning 
the BLUE units. 
18 Ship number 
J=1 Sum of the detection ranges on raiders 


J=2 Sum of the detection ranges on missiles 
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APPENDIX E 


GLOSSARY FOR PROGRAM VARIABLES 


a. MAIN PROGRAM 


ALFA REAL*8 variable used for conversion of degrees to 
radzvans: 
AMF small real constant used for conversion from float= 


ing point variables to integer variables. 


ASMT The table containing RED ASM/SSM launch and target- 


ane anformatiions 

BWJ Beam width of the jammer. 

BWR Beam width of the radar receiver. 

HISA Vector used to store data for use with subroutine 
HESTG: 

HISO Vector used to store data for use with subroutine 
HISTG. 

HISTG Subroutine in the NPS computer library that pro- 
vides a histogram display of an input vector. 

HISX Vector used to store data for use with subroutine 
HISTG. 

IBAT An array used to keep track of the number of mis-— 
Siles remaining in each missile battery of a BLUE 
unate 

IFAG Integer constant containing the number of RED mis- 


Siles engaged by the BLUE fleet. 


IFBG Integer constant containing the number of RED 
raiders engaged by the BLUE fleet. 


IHIT Vector containing the number of hits received by 
each BLUE unit. 


TIF The number of RED missiles that failed in flight. 
IKS The number of RED missiles shot down by BLUE SAMs. 
ILF The number of RED missiles that failed when launched. 






he 
Loh pa es 











alogy 
ok? dete La 


. “w9g 
Obl 


ja Yoetoe¥ 


77 ay) 


; i ‘to 7oe¥ 
, 488 


ji ad? 


ima odf 


it, 
- 






‘) eis ‘Qi t 


oo 


ILS The number of RED missiles that were in flight at 
- the end of the iteration. 


IMF The number of RED missiles destroyed when the 
launching aircraft or ship was destroyed. 


IMK The number of RED missiles destroyed by SAMs over 
all iterations. 


IMM The number of RED missiles destroyed prior to launch 
over all iterations. 


INSECT The original seed value used with the random number 
generator. Used to reproduce the random number 
Chaan af XC 54). is set toy 1. 


IRT The starting address in the X vector for the next 
five input quantities. 


ISAG An array containing the number of RED missiles 
first engaged by each BLUE ship. 


ISBG An array containing the number of RED raiders first 
engaged by each BLUE ship. 


ITS The number of RED missiles that impacted on their 


taraat 
a oo} we 


JACK The array containing the BLUE ship detection and 
engagement counters. 

KOUNT The iteration number. 

fee Lhe, vector containing the data,for the hit. fre- 


quency histogram. 


NADAR The array that stores the BLUE radar and missile 
system availability for each iteration. 


NASM The number of ASM/SSMs. 

ND Index used to convert true bearings to cartesian 
bearings. 

NI The number of iterations to be performed. 

NMIF The number of missiles in the BLUE fleet. 

NR The number of armed and Giemence RED raiders. 

NRTGT The vector containing the RED units target list 
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NS 
NSAM 
PLIST 


RAID 


RAT 


RXATA 


SAMASM 


SETIMX 


SHIP 


SMAT 


STATS 


ae 


The number of BLUE ships. 
The number of BLUE missile salvos fired. 
The array containing the target lists. 


The table containing the RED raider position and 
weapon information. 


REAL*8 variable used to convert degrees to radians. 


Dummy array used to read in data and fill the X 
VECtor. 


The array containing the BLUE missile RED target 
engagement information. 


The subroutine that sets the internal clock if 
mC Soo! tS §Set. to. a. 


The table containing the BLUE ship position infor- 
mation. 


The table containing the BLUE missile battery en- 
gagement data. 


2. SUBROUTINE EFFECT 


IUNIT 


NEXTEV 


TIME 


TIMES 


The code containing the units involved in playing 
the next event. 


The next event type. 
The time the next event is to be played. 


The mternal clock “time wsed*for entry and: exit 
a Coo) 1S set. 


3. SUBROUTINE JULIE 


BNG 


CHANEL 


IA 


A REAL*8 variable used to convert degrees to radians. 


The array that keeps track of which BLUE missile 
channels are busy. 


tines xX vector address of the fire control radar 
class number. 
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IBAT 


IC 


IFAG 


IFBG 


PPIRST 


IFLASK 


IFOG 


IGAGE 


EIT 


IJ 
IOR 
Ip 
Is 


ISAG 


ISBG 


LT 


ITEMS 


IVENTS 


JACK 


The array that keeps track of the number of BLUE 
missiles remaining in each battery. 

ihe weadarz class iof, the fire control madar. associated 
with the number one director of the indicated BLUE 
missile battery. 


Integer constant containing the number of RED mis- 
Siles engaged by the BLUE fleet. 


Integer constant containing the number of RED 
raiders engaged by the BLUE fleet. 


Mia used ior iheckine, a, target! s priority. 
Fae used) to set the time by which all RED raiders 
must. reach Gheiy launch point, in order to launch 


their missiles. 


Flag to indicate a) horizon search mode, and 
b) the system has a launcher that is operable. 


A vector that keeps track of which targets have 
been engaged by BLUE missiles. 


A vector that keeps track of how many hits each 
BLUE unit has taken. 


ine whip lass Starting index in the X vector. 
ihe. Lard class .Sstartines index <n. the, X vector. 
Missile system class number. 

The ship class starting index in the X vector. 


Integer constant to determine how many RED missiles 
the ship was the first to engage. 


Integer constant to determine how many RED raiders 
the ship was the first to engage. 


Search radar class number. 


The number of events that have been retrieved from 
the events list. 


The array containing the events list. 


The array containing the ship detection and engage- 
ment statistics. 
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JAM 


KFOG 
LOVER 


MAMBAT 


MISLFT 


NADAR 


NASM 


NONO 


NR 
NS 
NSAM 


PROB 


RAID 


RELOAD 


RGE 


SHIP 


SMAT 


STATS 


TGTAV 


The vector that keeps track of the number of jammers 
operating on a given radar frequency band. 


Flag to indicate that a search radar is available. 
The number of events in the events list. 


The number of missile batteries available to the 
BLUE force. 


The array that keeps track of how many missiles a 
RED raider has left on board. 


The array that keeps track of the operational mode 
of the BLUE equipment. 


The number of RED missiles. 


The array that keeps track of engageability of a 
RED unit or missile with respect to each BLUE unit. 


The number of armed and unarmed RED raiders. 
The number of BLUE units. 
The number of BLUE missile salvos. 


The probability that a given piece of equipment is 
operational. 


The table with the raider location and weapon in- 
formation. 


The array that stores the time the launcher will 
bere loaded sby "ai tera missile firing. 


The launch range of the longest ranged missile 
carried by each raider. 


The table containing the ship position information. 


The table with the BLUE missile system launch in- 
formation. 


The table of BLUE detection ranges. 


The table that tells when a target will be avail- 
able for engagement by another BLUE unit. 
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4. SUBROUTINE EVENT 


ICODE The participating units-event tyne code that is 
Svored oan the events: fast. 

ID The event type. 

iT The time of the event converted to the integer 


Lormac, USed —b0 store 1t in the events List. 


ITB Used to determine where in the events list the next 
event to be stored should go. 


ITEMS the number of events retrieved. 

ITIM ihe time: as istered in, the events list. 

ITM The time as stored in the events list. 

IUNIT The units involved in playing the next event. 


IVENTS The events list. 


IVET The participating units as stored in the events 
List 
LOVER The number of events stored in the events list. 


NEXTEV The next event type. 
NUMBER ihewcode for the partireipartine units. 
i The time of the event. 


TIME The time of the event. 


5. SUBROUTINE REMOVE 


IUNYT The code stored in the events list of the event to 
be removed from the events list. 


IVENTS he events list. 


LOVER The number of items in the events list. 


6. SUBROUTINE TARGET 


IFIRST The flag that indicates that the counter has been 
zeroed. 


rT The game time in events list code. 
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IUNIT 
IVENTS 


MORE 


NPGS 


TYME 


The code of the event to be compared by PRIOR. 
The events list. 


The flag indicating that there is an event with 
the same playing time and of type three. 


The counter that is used as an index in searching 
the events list. 


Present game time. 


i. SUBROUTINE SUBEL 


ALT 

AREA 
DELR 
FCRH 


FRG 


HFCSR 


HORZ 


HRAG 


HRFC 


HSFR 


HT 


IA 


IC 
ICH 


ICL 


The raider's altitude. 


The fourth root: of the raider's radar cross section. 


The detection delay. 

tie cagar nOrttzon.or the tire control radar. 
The minimum of the clear environment detection 
range, the radar horizon, and the radar's theo- 
retical maximum detection range. 


Lae radar horizon of the fire eontrol radar. 


The combined radar horizons of the radar and the 
target. 


The combined radar horizons of the radar and the 
target. 


The combined radar horizons of the radar and the 
target. 


The minimum of the clear environment detection 
range, the radar horizons, and the maximum theo- 
retical detection range. 


The. height of the fire control radar antenna. 


Used to check the class numbers of the various 
systems. 


Missile system class number. 
Used to enter the NADAR table. 


Fire control radar class number. 
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ICPA 


IFA 
IFC 
IFCR 
IF 
IR 
ISI 
JACK 


JAM 


JAMM 


JCPA 


JFOG 


KEY 


KFLAG 


LCPA 


MFOG 


MJL 
MTY 
NADAR 
NASA 
NASM 
NB 


ND 


Flag telling whether the target will come within 
the maximum detection range. 


Integer to store the radar number. 

Integer to store the battery number. 

Fire ‘control radar class number. 

Radar frequency band. 

The raid number. 

The ship class starting index in the X vector. 
The table of engagements and detections. 


The table that keeps track of how many jammers 
there are on each radar frequency band. 


The subroutine for determining the detection range 
of a monostatic radar in a jamming environment. 


ihe five that tells whether or not the target will 
come within the detection range. 


The flag that tells whether or not this is the 
second bistatic detection zone. 


The type of probability distribution associated 
with this unit's time delays. 


Flag that says that this unit has a search radar 
that is operational. 


Flag that tells whether or not a target will close 
to within the detection range. 


Counter giving the number of second bistatic de- 
tection zones to be processed. 


Index for table lookup. 

Missile system class number. 

The table containing the equipment availability. 
The flag indicating that sadder is an aircraft. 
The number of ASM/SSM launched or destroyed. 

The battery number. 


The number of bistatic detection zones. 
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NF The director number. 


NFOG Flag telling whether or not there are any more 
bistatic detection zones to check. 


NJ The number of jamming raiders. 

NONO The table that contains the reason for non- 
engageability. 

NSR The search radar index. 

NSURF Flag that says that this ship has a surface to 
surface missile system. 

PROB Probabitity..of acgquis»tien. 

RAGJ Detection range in a jamming environment. 

RAID The table with raider location and weapon informa- 
tion). 

RCPS Range at CPA squared. 

REACT Subroutine for determining reaction time delays. 

RECJ Fire control radar detection range in the jamming 


environment. 


RG The minimum of the clear environment detection 
range, the radar horizon range, .aand the theoreti- 
cal maximum detection range. 


RGE Detection range. 

RGF Jamming environment detection range. 

RGFC Detection range of the fire control radar. 

RGI Detection range in a jamming environment. 

RGIA Storage for second bistatic detection range. 
RGIB Storage for second bistatic detection range. 
RGIC Storage for second bistatic detection range. 
RGM Missile range. 

RIP The maximum detection range, or present target 


position, whichever is smaller. 


RRFC Clear environment fire control radar detection 
range. 
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RRH 


R2 


SALT 
SAMVEL 
SHIP 
SIGMA 
SRAG 


SRFC 


STATS 


SVEL 


THR 
THSCAN 
TI 
TIME 
TIML 
TIMLA 
TIMLB 


TIMLC 


Radar horizon of the search radar. 


The detection range in the second bistatic detec- 


tion Zone. 

The radar horizon of the raider. 

The BLUE missile's speed. 

The, tabWeiwith® the: Shap position information. 
The raider radar cross Section. 

Missile range. 


Fire control radar clear environment detection 
range. 


Ship's detection range storage. 

Missile speed. 

Time raider will close to the given range. 
The designation time delay. 


Mhoa taimna of detect 
aA 44 v ws Maeve vl wv 


asta 


The time of detection in the second bistatic de- 
tection zone. 


The last time to fire. 
The firing time delay. 


imme the tarcet, wall close to fire control. 
acquisition range. 


Theoretical radar maximum detection range. 
Horizon scan detection delay. 

Time the target will close to the given range. 
Present game time. 

Time the target will go beyond the given range. 
Storage for lose detection time. 

Storage for lose detection time. 


Storage for lose detection time. 
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TIML2 Time of lose detection in the second bistatic de- 
; tection zone. 

TLS Last time to fire. 

TRGE Time of detection. 

TSHOOT Time of launch. 

TSR Total reaction delay. 

aA Storage for detection time in second detection zone. 

LTS Storage for detection time in second detection zone. 

TC Storage for detection time in Second detection zone. 

TWDS Reaction time delay for weapons designation equip- 
ment. 

VEL Raider speed. 

XCP The X coordinate of the CPA. 

HCP The Y coordinate of the CPA. 

8. SUBROUTINE SUBE2 

ALT Altitude of the RED ASM/SSM during flight in feet. 

ALTN2 Altitude of the RED ASM/SSM squared in miles. 

AREA The fourth root of the radar cross sectional area 
of the RED missile. 

ASMT The table containing RED missile launch and target- 
Ing i iniormat ion. 

DELR The detection time delay. 

DELTAX The distance in the X direction between the target 
ship and the ASM/SSM launch point. 

DELTAY The distance in the Y direction between the target 
Ship and the ASM/SSM launch point. 

DIST The three dimensional distance from the ASM/SSM 
iguneh=point “EoeGhe target ship. This distance 
assumes launch at cruising altitude. 

HORZ The combined radar horizons of the firing ship and 


the RED missile. 
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HRAG The combined radar horizons of. the firing ship and 
the RED missile. 


HT The BLUE radar antenna height. 

TA Fire control radar class number. 

LC The BLUE missile system class number. 

ICL Fire control radar class number. 

ICOP Constant for X vector entry set by SUBE2. 

ICPA Flag telling whether or not: a RED missile will come 


within the detection range. 


IDP Index counter used to make up the RED raider's 
temporary target list. 


IFA The director number. 

LEC The battery number. 

IFCR BLUE missile system class number. 

IFLAG aa madiecatang  REDt unas finthnes two missile 
Sau a@)e 


IFLASK Flag indicating the first air launch has occurred, 


set TMAX. 
IHIT Table of hits taken by BLUE shee. 
IJ Raider class number. 
IJI Radar frequency band. 
IOU Target ship number. 
IPQ Target hit number. 
IR Raider number. 
ISI Shap class starting index in the X vector. 
JACK Ship detection and engagement table. 
JAM Table that keeps track of the number of jammers on 


the radar frequency bands. 


J AMM Subroutine to determine monostatic radar detection 
ranges in a jamming environment. 
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LNUM 


MARIE 


MFOG 


MISLFT 


NADAR 
NASM 
ND 
NEXT 


NFOG 


NJ 
NMIS 


NONE 


NONO 


NOW 


Flag indicating that there are. more second bistatic’ 


detection zones to evaluate. 


Flag indicating that there were two bistatic detec- 
tion zones. 


Flag indicating that a second type ASM/SSM is being 
processed. 


Type of probability distribution to be used in 
computing the reaction time delay. 


Flag indicating BLUE unit has an operational search 
radar. 


Flag for one or two missile salvos. 


Piacasndacatine whether or not raid could target a 
particular BLUE unit. 


The class number of the raider's long range missiles. 


Number of hits required to put the BLUE unit out of 
2Ct1ON. 


Counter indicating the number of second bistatic 
detection zones, 


Array to indicate the number of ASM/SSM left on 
board ta trader. 


The array that stores BLUE Sadapment operability. 
The number of ASM/SSMs processed. 

The number of bistatic detection zones. 

Flag indicating the short range RED missiles. 


Flag indicating there are more detection zones to 
process. 


The number of jammers. 
The number of ASM/SSM of this type to be launched. 


The) index used) to draw random numbers from the 
random number vector. 


The array specifying the target-firing ship rela- 
casos hap 


The type missile being processed now. 
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NR 
NRTGT 


NS 


The number of raiders. 


The target list assignments. 


The number of BLUE units. 
The total number of ASM/SSM on board. 
The target lists. 


the probability of acqurszrtion in horizon scan 
mode of operation. 


The vector of (0,1) random numbers used to deter- 
mine missile launch times. 


The search radar detection range in a jamming en- 
vironment. 


The table with the raider location and weapon in- 
formation. 


The slant range to CPA. 


Phe RED missles anilight’ relvability.: 


The minimum of the clear, jamming, and theoretical 
detection ranges. 


Detection range. 

mecee ion range. 

Fire control radar detection range. 

Detection range. 

Storage for secondary bistatic detection range. 
Storage for secondary bistatic detection range. 
Storage for secondary bistatic detection range. 
BLUE missile range. 

Detection range squared. 

Minimum of present range and detection range. 


Horizon search clear environment detection range. 
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RRH 


RT 


R2 
SALT 
SAMVEL 
SHIP 
SIGMA 
SPD 
SRAG 


SRFC 


SRGE 


STATS 


SVEL 
TAC 
TD 
TDES 
TDET 


TD2 


TFIR 


TFIRE 


TFR 


THR 


THSCAN 


Radar horizon. 


‘Cumulative probability of targeting this target or 


one previous. 

Second bistatic detection range. 

The radar horizon of the RED weapon. 

BLUE missile speed. 

The table containing ship position information. 
The RED weapon radar cross sectional area. 

The RED weapon speed. 

Search radar detection range in a clear environment. 


Fire control radar detection range in a clear en- 
vironment. 


BLUE ship missile range. 


Array containing the sum of a ship's detection 
ranges. 


BLUE missile speed. 
Acquisition time delay. 
Time of detection. 
Designation time delay. 
Time of detection. 


Time of detection in the second bistatic detection 
zone. 


Last time to fire based on search radar range and 
maximum range missile. 


Firing time delay. 


Time target will come within the fire control de- 
tection range. 


Maximum theoretical radar range. 


Reaction time delay for horizon scan mode of opera- 
tzToeny 
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TIMF 
TIML 
TIMLA 
TIMLB 
TIMLC 
TIML2 
TINT 


TLS 
TM 
TMAX 
TRGE 
TSHOOT 
TSR 
TTA 
TTARG 
TTB 
TTC 
TWDS 
TX 


VEL 


Time target will come within the search radar de- 
tection range. 


Time target will come within the horizon scan 
detection range. 


Present game time. 


Time target will exit the fire control radar detec- 
tion range. 


Time the target will go beyond the search radar 
detection range. 


Storage for second bistatic detection zone loss of 
target time. 


Storage for second bistatic detection zone loss of 
target time. 


Storage for second bistatic detection zone loss of 
target time. 


Time target will go beyond detection range in the 
second bistatic detection zone. 


The vw timnel antervall overs which alll missitessimust be 
launched. 


Last tame to fire based on fire control radar range. 
Proposed RED missile launch time. 

Time. limit on launch of RED missiles. 

a) Hazard zone range squared, b) time of-detection. 
Time of proposed BLUE missile launch. 

Total reaction time delay. 

Storage for second detection zone detection time. 
The temporary target list. 

Storage for second detection zone detection time. 
Storage for second detection zone detection time. 
Delay for weapons direction equipment processing. 
Time between RED missile launchings. 


RED missile speed. 
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XCPA 
WEPA 


Zi: 


CPA position in the X direction. 
CRA wWosautionr in, the Y direction. 


(0,1) random number. 


o. SUBROUTINE SUBES3 


ASMT 


CHANEL 


EBTJR 
FT 


IBAT 


“COT 


NCH 
NC1 


NMS 


The table containing the RED missile launch and 
targeting information. 


The array that indicates which system channels are 
in use. 


The Y coordinate of CPA. 
BLUE missile flight time. 


Array giving the number of BLUE missiles available 
by battery. 


Flag indicating channel assigned. 
Ship class starting index in the X vector. 


Input code giving target, shin, svstem, and direc- 
tor antormation. 


A channel number. 

RED missile number for table lookup. . 
Firing ship number. 

Target ship number. 

Number of original target. 

Flag used to check director scheduling. 
Target number after searching PRIOR. 
Array to check launcher operability. 
Director number. 

Channel number. 

Number of channels associated with this radar. 
Channel number. 


Missile system number. 
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NSAM 
PROB 

| RAID 
REASON 
RELOAD 
RGE 
RJL 
RLT 


SAMASM 


SHIP 
SMAT 


TCF 


BLUE salvo number. 

Probability of kill, 

Table of raider positions. 

Code for non-engageability. 

Table that tells if the launcher is loaded. 
Systems maximum range. 

The X coordinate of CPA. 

Launcher reload time. 


The table of BLUE missile and RED target engagement 
data. 


Table of ship positions. 
Table of BLUE missile system data. 


Dummy vector for storing the times the channels 
will be free for other use. 


Time director first available for channel use. 
Time director needed. 


Time target available for engagement by another 
Una t. : 


Time channel free. 

Present game time. 

Time director no longer needed. 
Time director needed for homing. 
Present target X coordinate. 


Present target Y coordinate. 


10. SUBROUTINE SUBE4 


fHIT 


Used as iaeclass-indicatorijecife: ASM type, or as an 
integer Constant, i.e. number of hits to sink the 
taneetaship: 


The table containing the number of hits sustained 
thus far in the game. 
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The input code containing the RED missile number 
and the BLUE target ship number. 


The target ship number. 

The RED ASM/SSM number for table lookup. 

The RED ASM/SSM number in the input code. 

The probabitsaty of Vsucecessful flight for the RED 


missile. Also used in determining the input and 
output time: of \ the subroutine’ if °K(55) is set 


11. SUBROUTINE SUBES 


ASMT 


CHANEL 


ID 


IFAG 


IFBG 


IGAGE 


TP 


IQ 
ISAG 


ISBG 


JACK 


JAM 


The table containing the data on the RED missile 
engagement. 


The table that keeps track of which BLUE missile 
channels are being used. 


An integer constant used for the starting address 
imine X=array ofthe BLUE firing ship's class 
inputs. Also used as the number of RED ASM/SSMs 
launched up to the present game time. 


Integer constant that totals up.the number of RED 
missiles engaged by the BLUE fleet. 


Integer constant that totals up the number of RED 
raiders engaged by tne BLUE fleet. 


Integer constant that keeps track of which targets 
have been engaged. 


Radar class’ of the BLUE fire cControl! radar’ asso= 
ciated with the engagement. 


BLUE missile system class number. 


Integer constant that totals up the number of RED 
missiles first engaged by this firing ship. 


Integer constant that totals up the number of RED 
raiders first engaged by this firing ship. 


The table that keeps track of the number of missiles 
and raiders engaged. by the firing ship. 


The table that totals up the number of jammers 
operating on a particular radar frequency band. 
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JFK Integer used in table lookup. 

K The BLUE firing ship. 

KEY Acmused inibhe Statement +=] CALL REACT (KEnY,4. TA) == 
where 4 is the type of delay, KEY is the class of 
the particular distribution within the type, and 
TA is the delay. 

M The BLUE salvo number. 


MISLFT The table that keeps track of the number of mis- 
Siles remaining on board the RED raider. 


N The number of missiles in the BLUE salvo. 

NA The RED ASM/SSM number for table lookup. 

NASM The number of RED ASM/SSMs processed thus far. 

NB The BLUE ship's director involved in the engagement. 
NC The channel number of the fire control system en- 


gaged in the intercept. 


NDU The starting address of the class associated with 
the raider involved in the intercept 


an 
ais 


NMS The number of the ship's missile system involved 
in the intercept. 


NT The target number, 1-90 for a raider, 91-190 for a 
RED missile. 


NTOT The number of missiles left on board the raider. 


PB The probability that the BLUE missile will receive 
the proper support during flacht. 


PROB The overall probability of a successful intercept. 


SAMASM The table that keeps BLUE salvo information. 


SMAT The table that keeps BLUE system information. 
TA The time delay for target acquisition. 

TF The time delay for firing. 

TK The time delay for kill assessment. 
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12. SUBROUTINE SUBEG 


CHANEL 


The table that keeps track of which BLUE missile 
channels are being used. 


ihe fLirane Ship. 


The input code containing the ship, direc-— 


tor. 


system, 
and channel number to be released. 


The director number. 
The channel number. 
The system number. 


The input and output time for the subroutine, used 
it. Mi So) aS=Seti. 


13. SUBROUTINE TIMMY 


ALT 


ASMT 


ASMTIM 


BETA 


CS 
DELTA 


DELTAX 


DELTAY 


DIS 


The cross range factor in the intercept contour. 


The angle between the target flight path and the 
iineroLesagcht tot the’ target. 


The target aititude. 


The table containing the data on the RED missile 
engagement. 


Subroutine to determine RED missile position at a 
given time. 


The downrange factor in the intercept contour. 


The angle between the BLUE missile flight path and 
the target's direction of movement. 


The target course. 

The reciprocal of BETA. 

Used for a) the X distance from the target ship to 
the RED missile, b) the present distance in the X 
direction from the firing ship to the target, and 
c) the distance in the X-direction from the firing 
Ship. tothe: intercept point. 


Same as DELTAX but in the Y direction. 


Range to the intercept point. 
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DLT Range to the intercept point. 


DX Distance’ in, the. X.direction: from) the, initials target 
position to the intercept point. Used if inter- 
Gept point is over the firing ship. 


DY Distance in the Y direction from the initial target 
position to the intercept point. Used if the in- 
tercepte point as directly over the firing Ship. 


EX The exponent that defines the super ellipse of the 
intercept contour. 

FT BLUE missile flight time. 

FZL Left hand limit of the firing zone. 

FZL1 Lett hand limit—of the second firing zone if a 


second exists. 
FZR Rrehc hand limit of the firing zone. 
FZR1 Right hand limit of the second firing zone. 


GAMMA Angle between the BLUE missile flight path and the 
Lime, Of ~Suoht ptosthe ttarsetlak daunch: 


GSINE The sine of GAMMA. 

ID RED missile type. 

IND Flag indicating the f epsigoaalinty 10d the intercept. 
IR RED missile's launching raid Gere ae. 

JOAN RED missile's target ship number. 


JULIA RED missile number for table lookup. 


K Firing ship number. 
KL Starting index for the BLUE missile class. 
KR Raider class number. 


LOCTIM Subroutine for determining raider position at a 
desired time. 


M Target number. 
NMS BLUE missile system number. 
NS The, number. of ‘ships adn the BLUE force. 
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PHI The BLUE firing bearing. 


PHIi The angle of axis rotation. for intercept contour 
engagement evaluation. 


Pr 3.141592654. 

PIZ The hazard range squared. 

PK The single shot kill probability. 

PM The probability of successful launch. 
PROB PK times PM. 


RADIAN A function for converting angles in degrees to 
angles in radians. 


RAID The table with the raider position and weapon 
information. 


RAT A REAL*8 variable for use with RADIAN. 


REASON A code to indicate the reason for non--engagement. 


RMAX The maximum range of the BLUE missile. 

RMIN The minimum range of the BLUE missile. 

RUTH Target altitude in miles. 

RXATA Dummy vector use to determine the firing zone. 
SEEK Y intercept of the intercept contour seeker limi- 


tation approximation. 


SLOPE The slope of the line that approximates the seeker 
Look Limitation. 


SSPD Speed of the BLUE missile. 


THETA The angle, measured in the normal X,Y plane, of the 
line from the firing ship to the present target 


Positron. 
fT The reciprocal of THETA. 
TIME The time of the proposed missile launch. 
PRL The time of intercept. 
Tx The time the target will enter the maximum missile 


range for this system. 
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14. 
AREA 
COFF 
DELX 
DELY 
HJBW 
HRBW 
ITS 


JJTBR 


The time the target will depart the maximum missile 
range for this system. 


The speed of the target. 
The X coordinate of the intercept. 


The X intercept coordinate after axis rotation and 
Lxeansilation. 


The X intercept coordinate after axis translation. 
The X coordinate of the present target position. 

The end of game X coordinate for the target if it 
isea raider, or the target ship. coordinate if the 


BLUE missile's target is a RED ASM/SSM. 


The initial position of the target in the X direc- 
taLOn: 


Des coordinate, of the intercept. point. 


The Y intercept coordinate after axis rotation and 
translation. 


The Y intercept coordinate after axis translation. 
The Y coordinate of the present target position. 
The end of game Y coordinate for the target. 

the Y coordinate of the initial target position. 


A sorting subroutine. 


SUBROUTINE JAMM 


Fourth root of ttargetiis rnadarncross: section. 
Jammer's main beam/side-lobe coefficient. 

Hareetis! track «change in X direction minus 0.1 nm. 
Marret'ismtirack «hanse an Yudirechion! minus 0.2 nm. 
One-half jammer beam width in radians. 

One-half radar beam width in radians. 

Target ship number. 


Bearing from jammer to the jamming target. 
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JRBRG 
LJBW 


LTBRG 


NASM 
NR 
PJ 
RCOF 


RCOFF 


RCPS 
RJBRG 
RJBW 
RJHS 


RRJS 
RTBRG 
SUMJP 
TBRG 
TIMC 
TIMI 
TIMX 
TIMY 


VD 


Bearing from the jammer to the detecting radar. 
Left bearing of the jamming beam. 


Left bearing of the detecting radar's beam while 
pointed at the target. 


Number of the ASM. 
Number of raiders. 


Power of the jammer on the detecting radar's fre- 
quency. 


Detecting radar's side-lobe coefficient. 


Detecting radar's side-lobe coefficient for specific 
jammer-target-radar geometry. 


Range to track CPA squared. 

Detecting radar to jammer bearing. 
Right bearing of jammer beam. 

Detecting radar-jammer horizon squared. 


Range from the detecting radar to the jammer 
squared. 


Right bearing of the detecting radar's beam width 
while pointed at the target. 


Sum of the jamming power arriving at the detecting 
radar. 


Target bearing from the detecting radar. 
Time constant used to adjust midpoint time. Base 
on target's launch time or game time of target's 


Initial position. 


Time target is at its midpoint in detecting radar's 
detection envelope. 


Lime tarcet 1s at the X coordinate of its midpoint 
in the detecting radar's detection envelope. 


Timeotareeteis atthe Y coordinate of its\)midpoint 
in the detecting radar's detection envelope. 


Divisor used in obtaining the X and Y components 
of the target's velocity. 
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Vx 
VY. 
XC . 


XCOF 


XJ 


XJJ 


XJJT 


XJR 


XJT 
XO 
XR 


XRT 


X1 
YC 


Yd 


YJJ 


YJJT 


YJR 


eT 
YO 
yR 


¥RT 


X component of the target's velocity. 
Y component of the target's velocity. 
Ri coordinate of target's CPA to detecting radar. 


Product of the jammer and detecting radar coeffi- 
cients based on the jammer-target-radar geometry. 


X coordinate of the jammer when the target is at 
its midpoint. 


X coordinate of the jammer at game time 0. 


Difference in X coordinates of the jammer's posi- 
tion and the jammer's target position. 


Difference in X coordinates of the detecting radar's 
position and the jammer's position. 


X coordinate of the jammer's target. 
X coordinate of the target's terminal position. 


X coordinate of the detecting radar's position. 
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Differencesan Y coordinates of the dete 
position and the target's midpoint posi 


Bake) 
x coordinate of the target's initial position. 


Y coordinate of target's CPA to detecting radar. 


Y coordinate of the jammer when the target is at 
itssmidpoint. 


Y coordinate of the jammer at game time O. 


Difference in Y coordinates of the jammer's posi- 
tion and the jammer's target position. 


Difference in Y coordinates of the detecting radar's 
position and the jammer's position. 


Y coordinate of the jammer's target. 
¥Y coordinate of the target's terminal position. 
Y feoordinate (oft hetdetect ing radars position. 


Difference in Y coordinates of the detecting radar's 
position and the target's midpoint position. 
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Y coordinate of the target's initial position. 


Target altitude in nautical miles. 


SUBROUTINE BIDET 


Coetitecient used in Solution of bistatie detection 
and radar horizon equations. See Appendix A for 
bistatic equation, 

2A. 

A squared. 

Similar to A. See Appendix A. 

Similar to A. See Appendix A. 


Constant used in solution of receiver and trans-— 
mitter radar horizons equations. 


Constant of bistatic detection equation. See 
Appendix A. 


Bistatic radar detection range in clear environ- 
ment against one square meter target. 


Array used only in calling subroutine DTIMD. 


See definition in subroutine JAMM portion of this 
glossary. 


Range squared from track CPA to bistatic receiver 
position. 


Range squared from track CPA to bistatic transmit-— 
ter position. 


Similar to A. See Appendix A. 
Same definition as in subroutine CPA. 
Same definition as in subroutine CPA. 
Similar to A. See Appendix A. 


Constant used in the solution of bistatic receiver 
and transmitter radar horizon equations. 


Defined in subroutine JAMM. 


Receiver's radar horizon range squared. 
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RCOF 


RCOFF 


RRH 
RRJS 
RTBRG 


SUMJP 


Transmitter's radar horizon range squared. 
Defined in subroutine JAMM. 

Flag set to 1 for stationary targets. 

Ship number of bistatic transmitter. 


Array indicating the number of active jammers on 
each radar frequency band. 


Flag set to 1 when the algorithm has completed the 
jamming section of the detection phase. 


Defined in subroutine JAMM. 
Defined in subroutine JAMM. 
Defined in subroutine JAMM. 
Defined in subroutine JAMM. 
ASM number. 

Defined in subroutine JAMM. 


Track change in Y direction divided by track change 


ann dareet 10M. 
Defined in subroutine JAMM. 


Defined in subroutine JAMM. 


Square of the range from the bistatic receiver to 
Statronary otarcet"s) position. 


Defined in subroutine JAMM. 
Defined in subroutine JAMM. 
Defined in subroutine JAMM. 


Raduealy part (of ‘solution to quadratic equation for 
time above bistatic receiver's radar horizon. 


Range to receiver's radar horizon. 
Defined in subroutine JAMM. 
Defined in subroutine JAMM. 


Defined in subroutine JAMM. 
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An array containing the bistatic detection times. 


Defined in subroutine JAMM. 


Time constant based on target's launch time of 
missile. Otherwise start of game time. 


Array containing the real parts to the solution of 
the bistatic equation. 


Array containing the complex parts to the solution 
of the bistatic equation. 


Defined in subroutine JAMM. 


Time target 
Time target 


Time target 
range. 


Time target 
range. 


Time target 


Time target 


enters receiver's radar horizon range. 
departs receiver's radar horizon range. 


enters transmitter's radar horizon 
departs transmitter's radar horizon 


enters bistatic radar detection envelope. 


iS) Liest above both transmitter and 


receiver's radar horizon. 


Time target 
velope. 


Time. target 


departs bistatic radar detection en- 


is no longer above both the transmitter 


and receiver's radar horizon. 


Target's velocity in nautical miles per minute. 


Defined in subroutine JAMM. 


X component 


Y component 


of target velocity. 


of tarecet velocity. 


Defined in subroutine JAMM. 


Defined in subroutine JAMM. 


Square of the range from the bistatic transmitter 
GQ the stationary target's position. 


Defined in subroutine JAMM. 


Defined in subroutine JAMM. 
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XJJT Defined in subroutine JAMM. 


XJR Defined in subroutine JAMM. 

XO X coordinate of target's terminal position. 

XR x coordinate of bistatic radar receiver pesition. 

XRAD Similar to RRAD above. 

XRC X coordinate of targetitrackoCPA to’ bistatic 
receiver. 

XRT Defined in subroutine JAMM. 

XX X coordinate of the bistatic transmitter. 

XXC Mecoordinate of target track CPA to bistatic 
transmitter. 

el. x coordinate of target's initial position. 

YS Defined in subroutine JAMM. 

YJJ Defined in subroutine JAMM. 

YJJT Defined in subroutine JAMM, 

YJR Defined in subroutine JAMM. 

YJT Defined in subroutine JAMM. 

YO Wicoordinate of.target's terminal position. 

YR Y coordinate of bistatic receiver's position. 

ERC Me coordinate on tarcet, track CPA to, bistatic 
receiver. 

YRT Defined in subroutine JAMM. 

YX Y coordinate of the bistatic transmitter position. 

¥XC Micoonrdinate: of chargé track «CPA tto ‘bistatic ‘trans- 
mater. 

Y1 Yicoordinateiofwtarcet isicindtial position’. 

ZR Bistatic receiver antenna height in feet. 

ZRN Bistatic receiver antenna height in nautical miles. 

ZX Bistatic transmitter antenna height in feet. 
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ZXN Bistatic transmitter antenna height in nautical 


miles. 
Z1 Target altitude in nautical miles. 
Z12 Target altitude in nautical miles squared. 


16. SUBROUTINE CPA 


DELX First used as a test variable to determine if track 
passes through reference position. Second usage is 
to indicate track change in X direction. 


DELY First used as a test variable to determine if track 
passes through reference position. Second usage is 
to Gnadweate track change in Y direction. 


Q Track change in Y direction divided by track change 
my x direction. 


mC X coordinate of CPA position. 

XI xicourdanate tof anit ial track posit lon. 

xo xscoordanate of terminal track position. 
XP X coordinare of reference position. 

aC Miccoords nate ot (CPA position. 

¥T Y coorGgsnate of Initial track Hoenn ont 

YO Y coordinate of terminal track position. 
YP Y coordinate of reference position. 


17. SUBROUTINE RCHECK 


DELX The distance along the X axis from the raider 
position to the potential target position. 


DELY The distance along the Y axis from the raider 
position to the potential target position. 


DIS The square of the distance from the raider to the 
potential target. 


HT The height of the highest radar on the potential 
target ship. 


HTR the tradarchonszon «ton (HT: 
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RADAR 


RADAR2 


RR 


er 


The class number of the potential target ship. 


the raider number. 


The potential target ship number. 


Flag for go/nogo situation. 


The combined radar horizons of the raider and the 
potential target. 


RADAR squared. 


The maximum theoretical range for the raider's 


radar. 


RR squared. 


Radar horizon. tor the raider. 


The table containing the BLUE ship location infor- 


mation. 


The raider's present position along the X axis. 


The ship's present position along the X axis. 


The raider's present position along the Y axis. 


The ship's present position along the Y axis. 


18. SUBROUTINE TMDASM 


B1 
C 


DELTAX 
DELTAY 
DELX 


DELY 


Coefficient A in the quadratic equation. 


Cociiveient BSB ian the ‘quadratic equation’. 


Ban ¥ = mx + b 


Coefficient C in the quadratic equation. 


The distance in 
tion to the end 


The distance in 
tion to the end 


The distance in 
position’ ‘to’ the 


The distance in 
position to the 


the X direction from present posi- 
of game position. 


the Y direction from present posi- 
of game position. 


the X direction from the intercept 
present position. 


the Y Girection from the intercept 
present position. 
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DEX 


DEY 


DIC 


DIF 


D1 


D2 


IFLAG 


IRAID 


ISAMBT 


The distance in the X direction from the firing 
ship to the target at present time, 


Hhe distance in the Y direction from the firing 
Ship to the target at present time. 


The distance in either the X or Y direction from 
intercept to end of game. 


The distance in either the X or Y direction from 
present position to end of game position. 


Distance from present position to first intercept 
point, squared. 


Distance from present position to second inter- 
cept point, squared. 


Flag indicating target is already within the indi- 
cated range. 


Raider number of target or launching raid. 
Battery number. 
Target ship number when target is an ASM/SSM. 


PIAS indicating taresetomovinge in a north or south 
darection. 


Raider number. 
Firing ship number. 
The raid array. 


The range to the target or to intercept, whichever 
is smaller. 


Bee-L ance 

hes ship array . 

min VYa=omk ib. 

Range of interest. 

Time target crosses range. inbound. 
Time target crosses range outbound. 
Present game time. 


Be = e4AC. 
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Target speed. 
Inbound crossing point. 


Inbound crossing point. 


19. SUBROUTINE TMDBOM 


A 
B 
BA 
C 


DELTAX 
DELTAY 
DELX 
DELY 
DIC 


Dir 


A in the quadratic equation. 
Bain the qiadratiic equation. 
ban Y= mx + b. 

Gan *the quadratic equation. 


The distance from present position to the end of 
game position in the X direction. 


The distance from the present position to the end 
of game position in the Y direction. 


Distance from the intercept point to present posi- 
ton an the Xx<{“dairection. 


Distance in the Y direction from the intercept 
Point CO Cue present position. 


Distance from the intercept point to the end of 
game position in the X direction. 


Distance from the present position to the end of 
game position in the X direction. 


Raider number. 

Rava array ; 

The launch range. 

mine Yr -mxX +b: 

Time raider will reach the launch range. 


The raider speed. 
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20. SUBROUTINE LOCTIM 


#D 


DELTAX 


DELTAY 


DL 


YZ 


The distance between the raider end of game posi- 
bron and the position stored in the RAID table. 


The distance along the X axis from the end of game 
position to the position stored in the RAID table. 


Distance along the Y axis. 

Distance traveled in the time interval defined as 
the difference between the time stored in the RAID 
table and the time for which the raider position is 
aesired. 

The raider number. 


The time stored in the RAID table. 


The table that centains the position and weapon 
information. 


The time for which the raider position is desired. 
The raider speed. 

the X position at time T. 

The X position stored in the RAID table. 

The Y position at time T. 


The Y position stored in the RAID table. 


21. SUBROUTINE ASMTIM 


ASMT 


DL 


MA 


The table containing the data on the RED missile 
engagement. 


The distance between the missile launch point and 
the impact point. 


the distance the missile will travel in the time 
defined as the difference between the launch time 
and the time for which the missile position is 
desired. 

The RED missile type. 


The RED missile target number. 


The RED missile number for table lookup. 
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time for which the ASM/SSM position is desired. 
missile speed. 
XK position at the tame TI. 


Y posiztaon «ate the: t amen T,; 
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